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Abstract 

In  this  work  we  study  the  pollution-error  in  the  h-version  of  the  finite  element 
method  and  its  eflfect  on  the  local  quality  of  a-posteriori  error  estimators.  We 
show  that  the  pollution-effect  in  an  interior  subdomain  depends  on  the  relationship 
between  the  mesh  inside  and  outside  the  subdomun  and  the  smoothness  of  the 
exact  solution.  We  also  demonstrate  that  it  is  possible  to  guarantee  the  quality  of 
local  error-estimators  in  the  interior  of  a  finite-element  mesh-patch  of  interest  by 
employing  meshes  which  are  sufficiently  refined  outside  the  patch. 


1  Introduction 

The  design  and  certification  procedures  for  aerospace  structures  require  an 
accurate  stress  analysis  capability.  Detailed  stress  analyses  of  complex  aircraft 
structures  and  their  subcomponents  are  required,  especially  during  the  certifica* 
tion  phase  of  the  aaal}rsis.  The  finite  element  method  has  become  the  method  of 
choice  for  the  analysis  of  complex  structures  such  as  aircraft,  nuclear  power  plant 
components  and  automobiles.  Before  manufacturing  a  component,  the  designer 
must  be  able  to  predict  its  behavior.  This  predictiim  is  based  on  a  formulation  of  a 
mathematical  model,  its  computational  analysis,  experiments,  and  experience  with 
existing  constructions  and  their  failures.  Because  of  various  uncertainties  which 
necessarily  occur,  the  goals  of  advanced  design  analyses  (in  aircraft,  nuclear  indus¬ 
tries,  etc.)  are  often  stipulated  in  the  design  codes  (which  are  rh^nging  over  time). 
The  question  of  the  principles  of  safety  is  directly  related  to  these  codes.  For  ex¬ 
ample,  in  the  design  code  USAF-MIL-A-83444,  used  in  military  aircraft  design,  it 
is  required  that  components  based  on  the  prindples  of  “non-inspectable  slow  crack 
growth”  must  be  designed  under  the  assumption  that, 

a)  the  as-fabricated  structure  contsdns  flaws  of  a  size  just  smaller  than  the  non¬ 
destructive  maximum  undetectable  flaw-size; 

b)  the  flaws  are  assumed  to  exist  in  the  form  of  crack-like  defects  with  most 
iinfavorable  location  and  orientation. 

The  design  code  requires  that  the  mathematical  formulation  and  its  computational 
analysis  must  reliably  and  conservatively  predict  both  the  sizes  of  the  growing 
cracks  and  the  residual  strength  of  the  component. 

As  a  result  of  the  requirements  of  the  design  code  many  times  it  becomes  neces¬ 
sary  to  analyze  a  small  portion  of  the  structure,  such  as  a  stress  critical  component, 
in  greater  detail,  e.g.  the  fuselage  of  an  aircraft  may  be  subjected  to  accurate  anal¬ 
ysis  in  the  neighborhood  of  structural  detaib  like  openings,  cracks,  etc.  A  standard 
procedure  is  to  analyze  the  entire  structure  using  a  global  discretization  which  is 
often  referred  to  as  global  analysis  and  to  further  analyze  the  structure  near  crit¬ 
ical  regions,  identified  by  the  global  analysis,  by  performing  local  analysts.  Such 
procedures,  which  are  known  as  global/load  analysis,  are  becoming  quite  popular 
in  engineering  practice. 

1.1  Generalities  about  global/local  analysis 

We  would  like  to  discuss  further  some  of  the  ideas  involved  in  the  technologies  of 
global/local  analysis.  In  the  literature  of  global/local  analysis  one  often  encounters 
the  terms  global  model  and  local  model.  We  would  now  like  to  comment  on  the 
meaning  of  the  word  model  which  is  often  used  with  different  meanings  in  the 
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engineering  community.  There  are  two  interpretations  of  the  finite  element  method 
(see  also  [1]).  The  first  one  considers  the  method  as  an  extension  of  the  standard 
method  of  structural  analysis  in  which  the  ori^nal  (physical)  problem  is  replaced  by 
another  (physical)  problem  which  is  finite  dimensional  and  can  be  computationally 
analyzed.  For  example,  in  the  late  50*8  it  was  common  practice  to  ‘‘model”  the  two- 
dimensional  elastic  continuum  by  a  system  of  trusses  connected  by  shear-panels. 
This  physical  “modelling”  method  was  called  finite  element  modelling.  The  second 
interpretation  is  to  understand  the  finite  element  method  as  an  approximation 
method  for  solving  mathematical  problems.  The  second  interpretation  is  the  main 
interpretation  today,  although  the  (misleading)  term  “finite  dement  modelling”  is 
stiU  used,  especially  in  engineering.  A  model  is  a  continuous  mathematical  problem 
and  is  approximated  by  a  finite  element  method  whidt  employs  meshes. 

It  follows  from  the  above  discussion  that,  if  the  same  mathematical  model  is 
solved  during  the  global  and  the  local  analysis,  the  terms  “global  model”  and  “local 
model”  are  misleading;  instead  the  terms  global  mesh  and  local  mesh  must  be  used. 
The  objective  of  global/local  analyses  of  this  type  is  to  obtain  accurate  i^proxima- 
tions  of  the  solution  quantities  of  engineering  importance  (stresses,  thermal  fluxes) 
in  local  regions  of  interest,  or,  in  other  words,  to  obtain  an  approximate  solution 
of  the  mathematical  model  with  error  less  than  a  specified  tolerance  in  the  region 
of  interest  (for  a  measure  of  error  vdiich  is  dictated  by  the  goals  of  the  analysis). 

It  is  also  possible  to  attach  another  meaning  to  the  term  global/local  analysis. 
For  example  a  3-D  (three-dimensional)  elastic  continuum  in  a  thin  domain  may  be 
analyzed  using  different  plate  models  in  different  regions  in  order  to  meet  desired 
levels  of  accuracy  with  respect  to  the  exact  solution  of  the  problem  of  3-D  elasticity. 
In  this  case  the  terms  “global  model”  and  “local  model”  may  be  appropriate.  This 
methodology  is  known  as  hierarchical  modelling  [1,  2,  3,  4,  5]  and  its  objectives 
are  to  approximate  the  solution  of  an  exact  mathematical  model  (for  example  the 
problem  of  3-D  elasticity)  by  employing  solutions  of  approximate  mathematical 
models  (e.g.  plate  or  shell  theories)  in  order  to  meet  specified  error  tolerances  with 
respect  to  the  solution  of  the  exact  mathematical  model,  in  the  regions  of  interest. 
The  approximate  modeb  are  continuoiis  mathematical  problems  which  are  derived 
frrom  the  exact  mathematical  model  using  certain  simplifying  assumptions. 

As  discussed  above,  the  term  global /local  analysis  has  been  employed  in  two 
different  contexts  in  practical  engineering  computations. 

(1)  In  one  context,  the  global  and  local  analyses  are  refined  or  enriched  mesh- 
techniques  often  employed  in  order  to  approximate  the  solution  of  a  single 
mathematical  model  throughout  the  domain.  An  example  of  this  type  of 
analysis  is  given  in  [6]:  In  that  work  a  methodology  is  outlined  in  which  a 
relatively  coarse  global  mesh  is  used  to  analyze  a  complex  structure.  The 
critical  regions  requiring  a  more  detailed  analysis  are  subsequently  identified 
and  interpolation  regions  are  then  defined  around  these  regions  and  an  inter¬ 
polation  procedure  is  employed  to  determine  the  boundary  conditions  for  the 
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local  meshes.  The  local  meshes  which  lie  within  the  interpolation  region  are 
independent  of  the  global  mesh  and  are  more  refined  than  the  global  mesh. 
In  some  analyses  of  this  type  the  discrete  problems  in  the  local  meshes  may 
be  also  fully  coupled  with  the  discrete  problem  in  the  global  mesh. 

(2)  In  another  context,  the  global  and  local  analyses  are  performed  by  solving 
different  mathematical  problems  in  different  parts  of  the  domain.  In  [7] 
such  an  approach  was  employed  in  order  to  resolve  crack-tip  fields  in  fibrous 
composites.  This  technique  which  in  [7]  is  called  the  mesh-superposition 
method  may  be  described  as  follows:  A  composite  plate,  which  may  include 
cracks,  is  first  analyzed  as  a  homogeneous  orthotropic  continuum  using  a 
finite  element  mesh  which  is  called  the  macro-mesh.  The  local  behavior  on 
the  scale  of  the  heterogeneous  constitumts  is  then  analyzed  using  a  separate 
finite-element  mesh  called  the  micro-mesh  which  employs  elements  which 
are  small  enough  to  reflect  the  micro  variations  in  material  behavior  and 
is  superimposed  on  the  macro-mesh  in  the  neighborhoods  of  the  crack-tips 
where  the  critical  behavior  is  expected. 

FVom  a  study  of  the  literature  it  is  dear  that  the  emphasis  of  the  research  in 
global/local  analysis  has  been  on: 

(i) .  Postprocessing  the  results  of  the  global  analysis  in  order  to  determine  the 

subregions  requiring  detailed  analysis; 

(ii) .  Generating  refined  and/or  enriched  meshes  in  the  local  areas  of  interest; 

(iii) .  Imposing  interface-conditions  at  the  interface  between  the  global  and  local 

meshes  (or  global  and  local  mathematical  modeb); 

(iv) .  Developing  solutions  strategies  for  the  linear  algebra  problems  which  result 

in  this  type  of  analyses. 

Very  little  attention  has  been  paid  on  ensuring  good  quality  (i.e.  accuracy)  of  the 
finite  element  solution  in  the  region  of  interest. 

In  this  paper  we  will  demonstrate  that  the  accuracy  of  the  approximate  solution 
in  a  subdomain  of  interest  de|>ends  on  the  relationship  between  the  mesh  used  inside 
and  outside  the  subdomain  where  the  detailed  analysis  is  required  or  in  other  words 
accuracy  of  the  local  analyses  depends  on  the  accuracy  of  the  global  analysis. 

To  illustrate  this  point  we  give  the  following  example:  Let  us  consider  Laplace’s 
equation  in  the  L-shaped  domain  shown  in  Figs,  la- Id  with  boundary  conditions 
consistent  with  the  singular  solution  u(r,9)  =  ri  sin(|)  in  the  infinite  wedge  (see 
also  Section  2).  We  assumed  that  we  are  interested  only  in  the  solution  in  the 
hollow  square  subdomain,  which  is  shown  with  thick  perigram  in  the  upper  right 
comer  of  the  domain  in  Figs,  la-ld.  We  computed  finite-element  solutions  of  the 
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Lapladaa  using  the  standard  displacement  method  and  we  employed  four  meshes  of 
linear  elements,  as  shown  in  Figs,  la-ld.  In  these  meshes  the  grid  in  the  subdomain 
was  kept  fixed  while  the  mesh  employed  outside  the  subdomain  was  coarsened  in 
order  to  illustrate  the  effect  of  the  outside  grid  on  the  acoiracy  of  the  finite^lement 
solution  in  the  subdomain.  In  Fig.  1.1a  we  show  an  adaptive  finite-element  mesh  of 
linear  elements  which  was  generated  using  a  remeshing  algorithm,  over  the  domain. 
For  this  mesh  an  accuracy  of  1.09%  for  the  relative  error  in  the  energy-norm  over 
the  subdomain  was  achieved  in  the  hollow  square  subdomain.  In  Figs.  1.1b,  1.1c, 
l.ld  we  show  grids  in  which  the  mesh  outside  the  subdomain  was  progressively 
coarsened  while  the  mesh  inside  the  subdomain  was  kept  fixed  as  in  Fig.  la.  Using 
the  grids  shown  in  Figs.  1.1b,  1.1c  and  l.ld  we  computed  approximate  solutions 
for  Laplace’s  equation  and  we  found  that  the  exact  relative  error  in  the  energy 
norm  over  the  subdomain  is  4.67%,  12.24%  and  23.9%,  respectively. 

FVom  this  example,  it  is  clear  that  the  accuracy  of  the  finite-element  solution 
in  any  subdomain  depends  on  the  mesh  outside  the  subdomain  and  that  in  order 
to  increase  the  accuracy  of  the  finite-element  solution  in  a  subdomain  of  interest 
one  has  to  employ  meshes  which  are  sufficiently  refined  inside  and  outside  the 
subdomain.  Hence,  in  general  we  cannot  say  that  a  refined  mesh  in  a  subdomain 
can  guarantee  a  better  solution  in  the  interior  of  the  subdomain.  The  situation  is 
very  similar,  but  more  complex,  in  the  case  of  hierarchical  modelling  and  will  not 
be  addressed  here. 

Based  on  these  initial  observations  it  is  clear  that  if  a  detailed  analysis  is  carried 
out  by  employing  a  refined  or  enriched  mesh  in  a  local  region  of  interest  with  no 
regard  to  the  mesh  outside  that  region,  no  gain  in  the  local  accuracy  may  be 
obtained.  Hence  the  idea  of  global/local  analysis  makes  sense  only  if  it  is  based  on 
algorithms  which  achieve  the  desired  accuracy  in  regions  of  interest  by  controlling 
the  relationship  between  the  inside  and  outside  meshes  in  order  to  employ  nearly 
minimal  total  number  of  degrees  of  freedom.  The  objective  of  this  work  is  to  make 
initial  progress  towards  this  goal. 

1.2  Introductory  remarks  about  error  estimation 

A-posteriori  error  estimation  is  the  only  rational  means  for  checking  the  local 
quality  of  approximate  solutions.  Therefore  it  is  important  to  know  under  what 
conditions  the  error  estimators  are  reliable  and  how  to  design  the  meshes  in  order 
to  guarantee  the  reliability  of  local  error-estimation. 

The  design  and  verification  of  error  estimators  for  approximate  solutions  of 
elliptic  boundary-value  problems  has  been  the  focus  of  intensive  research;  see  for 
example  [8-43]  and  the  citations  in  these  papers.  Studies  of  the  factors  which  affect 
the  local  quality  of  error  estimators  are  given  in  [41-45].  These  factors  are: 

1.  The  geometry  of  the  grid:  By  geometry  we  mean  the  topology  of  the  grid 


6 


7 


(b) 


Figure  1.1  (continued) 
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(c) 


Figure  1.1  (continued) 
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and  the  distortion  of  the  elements;  the  geometry  has  to  be  understood  in 
connection  with  the  differential  operator  (e.g.  [41]). 

£.  The  smoothness  of  the  solution:  The  smoothness  is  characterized  by  the 
regularity  of  the  solution  (see  [46]- [56])  which  depends  on  the  differential 
operator,  the  loads,  the  geometry  of  the  boundary  of  the  domain  and  the 
geometry  of  the  material  interfaces.  For  further  details  about  the  smoothness 
of  solutions  of  linear  elliptic  boundary-value  problems  in  polygonal  domains 
see  [49]-[56]. 

S.  The  pollution  error:  By  pollution  error  we  characterize  the  errors  in  the 
approximation  which  emanate  from  singiiiar  points,  boundary-conditions, 
mesh-  and  material-interfaces  and  which  may  affect  the  order  of  convergence 
of  the  approximation  throughout  the  domain  (see  [57]- [66]). 

In  [44]  and  [45]  we  presented  a  theoretical  framework  and  a  computer-based 
approach  for  checking  the  local  quality  of  estimators.  The  methodolc^  in  [44],  [45] 
cmables  one  to  obtain  a  quantitative  measure  of  the  quality  for  any  error  estimator 
for  interior  mesh-cells  uiq  in  the  interior  of  the  grid  (an  example  of  an  interior 
mesh-cell  is  given  in  Fig.  1.2;  the  mesh-cells  are  patches  from  the  mesh  which  may 
consist  of  a  few  elements,  possibly  one).  In  [45]  we  studied  the  robustness  of  several 
error  estimators  (which  are  used  in  practice)  for  interior  mesh-cells  for  the  complex 
grids  which  are  employed  in  practical  engineering  computations  and  we  identified 
several  robust  estimators. 

The  theoretical  framework  of  the  methodology  for  checking  the  local  quality 
of  estimators  (which  was  employed  in  the  study  of  estimators  in  [44]  and  [45]) 
assumes  that  the  pollution-error  in  the  mesh-cell  Uq  is  negligible  with  respect  to 
the  error  in  the  best-approximation  ddined  over  a  slightly  bigger  mesh-cell 
which  includes  Uq  and  a  few  mesh-layers  around  it.  This  assumption  implies  a 
proper  relationship  (which  depends  on  the  exact  solution  and  hence  on  the  data  of 
the  problem)  between  the  mesh-size  inside  and  outside  the  mesh-cell  Uq.  In  this 
work  we  show  that  this  assumption  is  necessary  (is  not  a  technicality  in  the  proofs 
given  in  [44])  i.e.  if  the  pollution-effect  in  is  large  (with  respect  to  the  error  in 
the  local  best-approximation),  nothing  can  be  said  about  the  reliability  of  any  error 
estimator  which  is  based  on  local  computations  in  Uq.  Moreover  we  demonstrate 
that  it  is  possible  to  control  the  pollution-error  in  any  interior  mesh-cell  by  proper 
design  of  the  mesh  and  hence  it  is  possible  to  ensure  the  reliability  of  the  error 
estimators  in  any  region  of  interest. 

In  this  work  we  will  study  the  effect  of  the  pollution  error,  due  to  singularities 
introduced  by  rough  domains  (domains  with  comer-points)  or  abrupt  changes  in 
the  boundary-conditions,  on  the  quality  of  error  estimators  for  interior  mesh-cells. 
Singularities  of  this  type  occur  in  typical  two-dimensional  geometrical  idealiza¬ 
tions  of  problems  in  solid  mechanics  and  heat-transfer.  Moreover  it  is  well  known 
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Figure  1.2.  (a)  A  finite-element  grid  with  an  interior  mesh-cell  (shown  shaded 
gray)  surrounded  by  the  mesh-cell  ^  and  three  mesh-layers  around  it;  (b) 
Close-up  view  of  the  mesh-ceUs. 
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(e.g.  [64],  [65])  that  singularities  of  this  type  pollute  the  approximation  throughout 
the  domain  while  singularities  introduced  by  rough  right-hand  sides  (concoitrated 
loads)  affect  the  approximation  only  in  the  neighborhood  of  the  singular  points. 
Therefore  the  pollution  error  caused  by  such  singularities  (due  to  rough  domains 
or  rough  boundary-conditions)  is  of  great  importance  in  practical  computations. 

In  this  work  we  will  analyze  the  pollution  error  due  to  rough  domains  and 
rough  boundary-data  in  the  h-version  of  the  finite  element  method,  its  effect  on 
local  error-estimation  for  interior  mesh-cells  and  how  the  pollution-error  can  be 
controlled  to  guarantee  the  reliability  of  local  error  estimation.  We  considered 
only  interior  mesh-cells  for  simplicity;  the  conclusions  of  this  study  are  also  true 
for  mesh-cells  which  abutt  to  the  boundary  of  the  domain  and  for  mesh-cells  near 
or  at  singular  points. 

The  main  points  of  the  work  are: 

1.  For  quasi-uniform  meshes  (often  used  in  practical  computations)  the  pollu¬ 
tion  error  could  be  significant  throughout  the  domain  depending  on  the  data 
of  the  problem  and  the  polynomial  degree  of  the  elements  employed  in  the 
mesh. 

2.  When  the  mesh  is  refined  locally  (in  a  region  of  interest  only)  the  pollution- 
effect  may  dominate  the  computation  and  there  may  be  no  gain  in  the  accu¬ 
racy  in  the  region  of  interest. 

3.  All  practical  error  estimators  are  based  on  local  computations  and  therefore 
cannot  detect  the  pollution  error. 

4.  When  the  mesh  is  nearly  equilibrated  in  the  energy-norm  the  pollution-effect 
in  the  derivatives  is  negligible. 

5.  The  pollution-effect  in  an  interior  subdomain  can  be  controlled  by  employ¬ 
ing  meshes  which  are  coarser  (outside  the  subdomain  of  interest)  than  the 
globally  nearly  equilibrated  grids. 

Following  this  Introduction,  in  Section  2  we  describe  the  model  elliptic  problem 
and  two  error  estimators  which  were  employed  in  the  numerical  examples.  In 
Section  3  we  give  numerical  examples  and  the  analysis  of  the  pollution-effect  for 
uniform  meshes  and  introduce  the  notions  of  near-field  and  far-field  pollution. 
In  Section  4  we  give  examples  that  demonstrate  that  when  the  grid  is  nearly 
equilibrated  in  the  energy-norm  the  pollution-effect  is  negligible.  In  Section  5  we 
demonstrate  that,  when  the  mesh  is  refined  locally  in  a  subdomain  of  interest,  there 
may  be  no  gain  in  the  accuracy  in  the  interior  of  the  subdomain  and  in  Section  6  we 
give  numerical  examples  which  indicate  that  it  is  possible  to  control  the  pollution- 
error  in  a  subdomain  of  interest  by  employing  grids  which  are  less  refined  (outside 
the  subdomain)  than  the  nearly-equilibrated  grids  (i.e.  the  globally-adaptive  grids). 
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2  Preliminaries 

2.1  The  model  problem 

Let  OCR*  denote  a  polygonal  domain  with  boundary  di)  —  F  —  U  f,-, 

where  F^  t  ss  1, . . . ,  are  the  open  straight  edges  connecting  the  endpcnnts  Af 
and  Ai^i  {Ai  —  Am+i)  FiS-  2*^)*  denote  the  internal  an^es  at  the 

vertices  by  (0  <  <  2x ;  if  s  2*  we  have  a  aUt-domain).  Let 

further  F  =  Fd  U  f  jv,  Fo  O  Fjv  =  0  where  Fd  is  the  Diriehlet-  and  Fs  is  the 
iVeumann-boundary. 

By  m  is  positive  integer,  we  denote  the  Sobolev  space  of  functions  with 

square  integrable  derivatives  up  to  order  m  with  the  norm,  seminorm 

ll-llir-co)  E  ll®‘“lln.  Nli-iO)  := 

0<k<m 

(®-«)(.)  ==  |(  E_l^«l’)(*)f .  == 

where tt :«  (01,03),  O(,  >0,  integers,  t  =  l,2,  |flr|  :s:oi  +  03.  Here  11*110 
denotes  the  L^-norm  over  (1. 

Let  T  {Ta}  be  a  ngvlar  family  of  triangulations  of  ft  (the  minimal  angle  of 
all  the  triangles  is  bounded  below  by  a  positive  constant,  the  same  for  all  meshes). 
The  meshes  are  not  assumed  to  be  quasi-uniform.  We  let  E  (resp.  Eint)  denote  the 
set  of  edges  (resp.  interior-edges)  in  the  triangulation.  We  will  consider  the  mixed 
boundary- value  problem  for  the  Lapladan: 


L(u)  :=  —V  •  Vu  ss  —  Att  =  0 

in  n 

(2.1.) 

u=  0 

on  Fn 

(2.14) 

du 

•—  :=  v«  'n=  g 
on 

on  Fn 

(2.1c) 

Here  u  is  the  solution;  n  denotes  the  extmor  unit-normal  on  Fjv ;  y  €  L^{Fn)  and 
is  analytic  in  each  Fj ,  j  s  1, . . . ,  Jlf . 

In  the  neighborhood  0{mf]  Tq)  =  fl  O  {«  |  |a;  -  Xg  |r3  <  ro  }  of  the  vertex  A^ 
{xi  is  the  position-vector  of  i4f),  the  exact  solution  can  be  writ^  in  the  form  (for 
homogeneous  boundary-conditions  in  the  ndghborhood  of  the  vertex  At) 
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(2.2.) 


tt(r/, 9t)  *  23 i^t) n*  +  Uo(r/, , 

with  imoother  than  the  fint  term,  where  are  the  polar  coordinatea 

with  leap^  to  the  vertex  Ag.  Hie  fhnctkms  and  the  exponents  aj  are 

{lin(ajtfi) ,  for  the  Dirichlet  or  mixed  problem  , 

(2.2h) 

coB(ajBe) «  f<v  the  Neumann  problem  , 


V>/  ’ 

(2i  -  l)x 

2^< 


for  the  Dirichlet  or  Neumann  boundary  omditions  , 


for  mixed  boundary-conditions  , 


(2.2c) 


(The  mixed  problem  here  means  homogeneous  Dirichlet  conditions  for  tfr  =  0  and 
homogeneous  Neumann  conditiims  for  $t  =  <pt)^ 

Remark  2.1.  hi  the  numerical  examples  we  considered  specially  unsmooth  sdutions 
of  the  type  given  by  (2.2a)-(2.2c)  because  such  solutions  occur  typically  as  solutions 
of  practical  boundary-value  problems. 

We  let  |u  €  |  u  ss  0  on  ri>|.  The  variational  formulatkm  oi 

the  model  problem  (2.1)  is:  Find  ti  €  such  that 

B(,u,v)=  (lUa) 

where 

Bo(».v):=  jf  Vii-Vt.,  »,  »€Jf.  (2.34) 

U  Tjy  =0  then  the  usual  restrictions  on  the  data  are  imposed  and  the  solution  u 
is  determined  up  to  an  arbitrary  constant.  We  let 

lll>’llls  ~  (Bs{v,v))i  =  (jf  |Vt.|»)*  (2.4) 

denote  the  energy-norm  over  S. 
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Let  ttfc  €  :=  jw  €  v  |r€  Vf{r)  V  r  €  T*,  v  =  0  mi  r^l  be  the 

piecewise  p-degree  finite  element  apprmdmation  of  the  solution  of  (2.1)  and  (2.2) 
which  satisfies: 

Find  Uj^  eSTr,  such  that 

^oK.v*)=  /  Vv*€5Jj.^  (2.5) 

•'»  nr 

The  error  is  u  *-  and  satisfies  the  rtaidwd-eqwition: 

Find  ej^  €  iffo  ^bat 

*>(«».»)  =  E  FAv)  V  V  e  (2.6) 

.er, 

where  Fr  :  Iff^,  R  is  the  residual  functional  given  by 

-PV(v)  :=  /  Wr  +  5  53  fvJ, ,  u  €  (2.7) 

Here  r,.  Au*  denotes  the  interior-residual  in  element  r  and  J(  is  the  inter¬ 
element-jump  in  the  normal  derivative 

lV(ufc  1^ )  -  V(u*  1^^,)]  •  *» .  e  € 

A:='2(s-^),  .CTw  (2.8) 

.0  ,  c  C  Fd 

where  n  is  the  unit-normal  assigned  to  the  edge  e  (in  an  arbitrary  but  unique  way) 
end  r^i ,  Ti.  are  defined  as  in  Fig.  2.2. 

2.2  The  error  estimators 

We  summarize  the  definitions  of  the  error  estimators  employed  in  the  numerical 
examples  below  (see  also  [44]  and  [45]). 

2.2.1  Implicit  element-residual  estimator 
Let  us  define  the  equilibrated  residual  functional 

Fr*(v)  -  Fr{v)  +  [  v9r,  V  €  HHt)  (2.9) 

JBt 
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Here  $r  €  (X*(dr))^  is  the  correction  of  the  edge-iesidunb  which  is  determined 
sttdi  thftt 


F^(v)  =  0  V  w  €  Vfir)  (2.10) 

Definiti<Mi8  oi  0^  which  satisfy  (2.10)  for  linear  elements  are  given  in  [14],  [15], 
[24],  [25],  [28];  in  the  numerical  implementatkms  for  linear  dements  we  employ  the 
definitkm  given  in  [14].  We  define  the  element  error-mdicator  function  Cr  as  the 
soltttkm  (which  is  unique  up  to  a  constant)  oi  the  local  problem: 

Find  e,  €  such  that 

»,)  -  (».)  V  €  B,  (2.U) 

Here  =  ^*(^)*  We  then  define  the  element  error-indicaton  ifr  for  the  energy* 
norm  of  the  error  by 

tlr  :=  \f^{Or,Or)  (2.12) 

In  the  computations  below  we  approximated  the  solutions  of  the  local  problems 
by  using  a  finite  dimensional  space  of  polynomials  of  degree  (p+ 1)  in  the  element, 
namely 


ST'  :=  {to  €  T>„..(t)  I  n;«-  -  0 }  (2.13) 

where  11^  :  (r)  is  an  interpolation-operator  defined  (see  [19],  [20]) 

over  the  element  r.  The  space  is  often  referr^  to  as  the  buhhle-spoee  (see 
[19]).  In  the  sequel,  we  will  call  £R-estimators  (element-residual-estimators)  the 
estimators  determined  from  (2.12)  by  computing  apprmdmate  solutions  of  (2.11) 
lErom  the  bubble-space. 

We  now  describe  how  the  corrections  for  the  edge-residuals  0^  are  constructed. 
We  let 


».i. = c  + ^(>0  •  '  £ 


where 


r-H,  if  r  =  Tta, 

Cr  =  {  (2.15) 

I  -1,  if  T  S=  T«., 

Here  it  is  assumed  that  an  edge-normal  n  has  been  assigned  to  the  edge  e  in  an 
arbitrary  but  unique  way  and  Ti.  and  are  defined  as  shown  in  Fig.  2.2.  Here 
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(2.16) 


where  A]^,  k=l,2  are  the  linear  8lu4>e*fimction8  defined  over  the  edge  e,  and 

Ai  .  *  =  1,2  ,  tCdr  (2.17) 

C 

The  integrals  in  (2.17)  are  obtained  from  the  linear  system  (see  Ladevexe 
and  Leguillon  [14]  and  Bank  and  Weiser  [15]  for  the  details) 

J  - -Fri^x)  ,  *  =  1 . (ne)x  (2.18) 

»r* 

where  9x  denotes  the  elementwise  affine  basis  function,  which  corresponds  to  the 
vertex  AT,  shown  in  Fig.  2.3b;  Tjf  denotes  the  ib-th  element  connected  to  the  vertex 
AT;  (ne)x  is  the  total  number  of  elements  connected  to  the  vertex  X. 

The  procedure  outlined  above  has  been  developed  by  Ladeveze  [14].  The  linear 
ssrstem  (2.18)  has  a  one  parameter  family  of  solutions.  Specific  choices  of  solutions 
are  suggested  in  [14]  and  [15].  In  the  numerical  implementations  we  employed  the 
choice  given  in  [14].  Below  we  ^ve  the  definition  of  the  edgewise-Unear  corrections 
6r  whi^  restilt  by  using  the  equilibration  procedures  of  [14]. 

Let  us  consider  the  vertex  X  and  let  Nx  denote  the  total  number  of  the  edges 
(or  the  elements)  connected  to  the  vertex  as  shown  in  Fig.  2.3a.  Let  us  also  denote 
by  Ti  and  e,-,  t  s  1, . . . ,  Nx  the  elements  and  the  edges  connected  to  X.  We 
determined  the  coefficients  ^**^1  which  are  associated  with  the  edge  e  and  the 
vertex  X  and  are  employed  in  (2.14);  here  the  indac<function  u{€,  X)  identifies  the 
node  X  in  the  local  enumeration  used  in  (2.14)  for  the  unknowns  associated  with 
the  edge.  The  values  of  the  coefficients  are  (see  also  [14]): 

Jt)  _  JT)  =  -L'£(N-i  +  l)F„  (*x)  (2.20a) 

^  ,=1 

(Ox)  ,  <  =  2,3,...,(Ar-l)  (2.204) 

The  edge-coefficients  are  computed  based  on  the  edge-normals  which 

rotate  counterclockwise  around  the  vertex  X  as  shown  in  Fig.  2.3a. 

Remark  B.2.  It  is  also  possible  to  define  the  equilibration  for  polynomial  spaces  up 
to  degree  p  (where  p  denotes  the  polynomial  degree  of  the  elements). 
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to  the  node  X\  (c)  The  local  enumeration  for  the  degrees  of  freedom  of  the 
correction  6  for  the  edges  connected  to  node  X . 


2.2.2  Error  estimator  based  on  patdi>averaging 

Here  we  outline  an  error-estimator  based  on  local  averaging  which  was  in¬ 
troduced  by  Zienldewicz  and  Zhu  [36,  37].  This  estimator  will  be  called  the  ZZ- 
estimator  in  the  examples  below. 

Let  uix  X—  (J  r'  denote  the  patch  of  elements  connected  to  vertex  X .  For 

X€ilf(r») 

each  patch  ux  we  recover  the  patch^projeetions  Wjf ,  by  solving  the  following  least- 
squares  problem: 


-  Vu,)  =  ^  4(<r  -  V».) ,  J(<,)  :=  ^  (».))’]  (2-21) 

Here  ,  m  =  1, . . .  ,nsp  denotes  a  set  of  sampling-points  in  the  patch  ux- 
The  element  error  indicators  for  the  ZZ-estimator  are  given  by 

I,,  :=  Ik”  -  Vu»||.  (2.22) 

Here  denotes  the  recovered  gradient  which  is  obtained  by  averaging  the  patch- 
projections  ax.  For  the  grids  of  triangular  elements  employed  in  this  study  a 
continuous  recovered  flux  a^^  is  obtained  by  taking  appropriate  averages  of  the 
flux  recovered  from  the  three  patches  associated  with  the  element.  In  this  work  we 
implemented  the  ZZ-estimator  for  elements  of  degree  p  —  1,  2,  and  3  as  follows: 

a.  Linear  elements  (p  =  1). 

The  nodal  values  of  a^^  are  recovered  from  the  patches  corresponding  to  each 
node.  The  recovered  C°-continuous  piecewise-linear  flux-field  a^^  over  the  domain 
is  constructed  as 


,”(*)  =  ♦,(,)  (2.23) 

Here  a\  is  the  recovered  flux  field  over  the  patch  ijtx  ^<1  value  of  the 

recovered  flux-field  at  the  vertex  X  and  is  the  piecewise-linear  basis  function 
associated  with  the  vertex  X. 

b.  Quadratic  elements  (p  =  2). 

A  piecewise  quadratic  Langrangian  representation  of  a^^  is  constructed  firom 

*=  ^*^x\x  ^  +  "“jraly)  ♦y(*)  (2-24) 

where  Y  is  used  to  denote  the  midside  nodes  for  the  edges;  ^ 

gran^an  quadratic  basis  functions  associated  with  the  nodes  Xj  Y  respectively; 
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Xiy  X2  are  the  vertices  at  the  endpoints  of  the  side  which  includes  the  midside 
node  Y. 

c.  Cubic  dements  (p  s  3). 

A  piecewise  cubic  Langrangian  representation  of  is  constructed  from 


■‘’I  \y  +  ®’Xa|y) 


X  "  *  y€{K,  .-y, ) 


(2.25) 


where  , . . . , are  used  to  denote  the  nodes  for  the  Lagrangian  cubic  element 
on  the  edges  and  Yf  denotes  the  interior  node  at  the  centroid  of  the  element; 
and  ^  cubic  TAtigrawgiAn  basis  functions  associated  with  the  nodes  X 
and  Y  respectively;  X,,  Ag  are  the  voices  of  the  triangular  element. 

Zienldewicz  and  Zhu  [36, 37]  proposed  two  types  of  least-square  problems  for  the 
recovery  of  namely  a  discrete  and  a  continuous  least-squares  patch-projection. 
Here  we  used  the  discrete  least-squares  patch  projection  and  we  employed  the 
sampling  points  shown  in  Fig.  2.4a,  Fig.  2.4b  and  Fig.  2.4c  for  linear,  quadratic 
and  cubic  triangles,  respectively. 


2.3  Definition  of  the  effectivity  index 

The  quality  of  an  error  estimator  in  a  mesh-cell  Hq  is  measured  by  the  effectivity 
index 


-J±. 


(2.26) 


Here  T|^  denotes  the  finite-element  mesh  which  includes  the  mesh-cell  in  its 
interior,  as  shown  in  Fig.  1.2,  ||efc||^  is  the  norm  (of  interest)  of  the  error  over 

(Oq,  is  an  error  estimator  for  this  norm  which  is  computed  in  terms  of  dement 
error-indicators  associated  to  every  element  r  of  the  mesh  T|^. 

In  [44],  [45]  it  was  demonstrated  that,  when  the  pollution  error  is  negligible,  the 
value  of  the  effectivity  index  in  the  mesh-cell  (Oq  is  influenced  (for  all  practical 
purposes)  by  the  geometry  of  the  mesh  in  a  mesh-patch  which  includes  and 
a  few  mesh-layers  around  it.  The  asymptotic  range  of  the  effectivity  index  in  the 
cell 
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Figure  2.4.  Location  of  the  sampling  points  employed  by  the  discrete  ZZ-estimator  in 
the  interior  mesh-patches,  (a)  Sampling  points  for  linear  elements;  these 
pmnts  are  located  at  the  centroids  of  the  elements;  (b)  Sampling  points  for 
quadratic  dements;  the  points  are  ocated  at  the  midpoints  of  the  sides  of 
the  elements;  (c)  Sampling  points  for  cubic  elements;  the  points  are  located 
at  the  vertices,  at  the  midpoints  of  the  sides  and  at  the  centroid  of  the 


elements. 


(c) 


Figure  2.4.  (continued) 
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(2.27) 


0<C^  <K^<C!^ 

is  defined  in  terms  of  the  constant  which  can  be  determined  by  a 

eomfuter^based  approach  (see  (44],  [45]).  Inequality  (2.27)  can  also  be  written 
in  the  form 

(2-M) 

which  expresses  the  equivalence  between  the  norm  of  the  error  and  the  estimator 
in  the  mesh-cell  u^. 

For  example,  for  linear  elements  (p  s  1)  and  for  the  Regular  mesh-pattem 

shown  in  Fig.  2.5a,  the  values  of  the  constants  are  (see  [44])  =  Cy  =■  1.00  for 

both  the  ER  and  the  ZZ  estimator  while  for  the  Criss-Cross  mesh-pattem  shown 
in  Fig.  2.5b,  the  values  are,  s  1.00,  1.414  for  the  ER  estimator  without 

equilibration,  s=  1.00  for  the  ZZ  estimator  and  the  ER  estimator  with 

equilibration.  FHirther  results  about  the  above  a'^d  other  estimators  for  various 
mesh-pattems  are  given  in  [44,  45].  From  the  :  results  it  is  clear  that  the  ER 
estimator  with  equilibrated  data  and  the  ZZ  are  robust  estimators  {(^  ns 

cif  »  1)  for  typical  meshes  used  in  engineering  computations.  As  it  was  pointed 
out  earlier  the  constants  C^,  can  be  determined  under  the  assumption  that 
the  pollution  error  is  negligible  in  a  neighborhood  of  the  mesh-cell  of  interest  (see 
[44]  for  the  mathematical  analysis).  In  this  work  we  demonstrate  that: 

1.  The  assumption  on  the  pollution  error  is  necessary  (i.e.  when  the  pollution 
error  is  significant  in  the  mesh-cell  (with  respect  to  the  error  in  the  best- 
approximation  in  u'^),  (2.28)  does  not  hold). 

2.  Robust  a-posteriori  error  estimation  in  a  mesh-cell  Uq  (and  hence  the  validity 
of  (2.28))  is  possible  by  controlling  the  pollution  error  in  Wq  through  proper 
design  of  the  mesh  outside  o;^. 
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(b) 


Figure  2.5.  (a)  Basic  mesh-oell  of  the  Regular  mesh-pattern  (b)  Basic  mesh-cell  firom 
the  Criss-Cross  mesh-pattem. 
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3  The  pollution-effect  for  uniform  meshes 

Practical  computations  in  elasticity  and  heat-transfer  are  often  performed 
using  meshes  which  are  nearly  uniform  (quasi-uniform).  It  is  therefore  important  to 
understand  the  pollution-effect  for  such  meshes.  Below  we  give  numerical  examples 
which  demonstrate  the  pollution-effect  for  uniform-meshes  of  linear,  quadratic  and 
cubic  triangles. 


3.1  Numerical  examples  of  the  pollution-effect  for  uni¬ 
form  meshes 

Here  we  present  numerical  examples  which  demonstrate  that  for  the  class 
of  uniform  meshes  and  the  class  of  specially  unsmooth  solutions  (with  algebraic 
singularities  of  the  type  given  in  (2.2a)-(2.2c))  the  pollution-effect  depends  on  the 
relationship  between  the  values  of  the  exponents  of  the  algebraic  singularities  and 
the  polynomial  degree  of  the  elements. 


S.1.1  Uniform  meshes  of  linear  elements 

We  considered  the  mixed  boundary-value  problems  for  the  Laplace  equation 
in  the  Irshaped  domain  shown  in  Fig.  3.1a.  We  applied  boundary- 

IT  1 

conditions  consistent  with  the  exact  solution  u(r,  0)  —  r®  sin(a^),  for  a  =  r —  =  - 
«  ,  2v>i  3 

^  2  /  1 

and  a  =  —  =  -  (fora=  -  homogeneous  mixed  boundary-conditions  (resp.  for 

u  V  o 

2 

a  —  homogeneous  Dirichlet  boundary-conditions)  were  applied  on  the  edges 
AiA2,  A^A^  which  emanate  from  the  singular-point  while  Neumann  boundary- 
conditions  were  applied  on  the  rest  of  the  boundary;  see  also  Section  2^ ,  and  we 

computed  finite  element  solutions  using  uniform  meshes  of  linear  triangles.  In 
Table  3.1  (resp.  Table  3.2)  we  give  the  effectivity  indices  for  the  mesh  layers 

n  =  1, . . .,  -  (a  tjrpkal  mesh-layer  is  shown  in  Fig.  3.1a)  for  ot  =  -  ^resp.  ®  • 

We  note  that  the  effectivity  indices  for  the  ER-estimator  and  for  the 
ZZ-estimator  are  practically  the  same  in  the  interior-layers.  We  also  note  that  for 

Ot  =  -  ^resp.  a  =  the  effectivity  indices  in  the  layers  at  a  fixed  distance  from 
the  singularity  converge  to  one  (resp.  zero)  as  the  mesh-size  k  tends  to  zero. 
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(«) 

Figure  3.1.  (a)  The  L-shaped  domain  with  boundary  -AjA, . . .  A^A^,  meshed  with  a 
nniforan  grid  of  triangular  elements  in  which  the  eighth  mesh-layet  is  shaded 
gray;  (b)  The  convex  domain  with  boundary  Aj  A, . . .  A4A1,  meshed  with  a 
uniform  grid  of  triangular  elements  in  which  the  eighth  mesh-layer  is  shaded 
gray;  (c)  The  triangular  domain  with  boundary  Aj  Aj . . .  A,Aj  medied  with 
a  a  uniform  grid  of  triangular  dements  in  which  the  sixth  mesh-layer  is 
shaded  gray. 
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Figure  3.1.  (continued) 
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Figure  3.1.  (continued) 
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Tkble  3.1.  (oomtinued) 
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l^ble  S.2.  (continued) 


Effectivity  index  for  endi  mesh*leyer 

Exact  solution  u(r 

,tf)  =  ri  Bin(|);  linear  elements 

Layer 

number 

n 

0.796 

0.678 

86.52 

0.853 

34.26 

0.666 

19.92 

0.572 

14.37 

0.519 

11.42 

6 

0.482 

9.57 

7 

0.453 

8.30 

8 

0.427 

0.429 

7.37 

9 

0.407 

0.409 

6.65 

10 

0.390 

0.392 

6.09 

11 

0.375 

0.375 

5.63 

12 

0.361 

0.362 

5.25 

13 

0.348 

0.349 

4.93 

14 

0.337 

0.337 

4.66 

15 

0.326 

0.327 

4.42 

16 

0.316 

0.317 

4.21 

17 

4.02 

18 

0.299 

3.86 

19 

0.292 

0.292 

3.71 

0.285 

0.285 

3.58 

21 

0.278 

0.278 

3.45 

22 

0.272 

0.272 

3.35 

23 

0.266 

0.266 

3.24 

24 

0.260 

0.261 

3.15 

25 

0.255 

0.255 

3.07 

26 

2.98 

27 

0.244 

0.244 

2.92 

28 

0.239 

0.239 

2.86 

29 

0.233 

0.234 

2.81 

0.228 

0.228 

2.76 

31 

0.221 

0.222 

2.74 

32 

0.215 

0.226 

2.72 
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S.1.2  Uniform  meshes  of  quadratic  elements 

We  OHuidered  the  mixed  boundary-value  problem  for  the  Laplace  equation  in 
the  ocmvex  dmnain  shown  in  Fig.  3.1b  and  ^plied  boundary-conditions 

jf  2 

consistent  with  the  exact  solution  u(r,fi)  s  r*sin(attf)  for  a  —  - —  =  -  and  a  = 

2(fij  3 

*■4/2 

—  —  X  (for  a  =  >  homogeneous  mixed  boundary  conditions  (resp.  for  a  = 

^  o  \  3 

4 

-  homogeneous  Dirichlet  boundary-conditions)  were  applied  on  the  edges 

Aj  A4  which  emanate  from  the  singular  point,  while  Neumann  boundary  conditions 
consistmt  with  the  exact  solution  were  applied  on  the  rest  of  the  boundary^  and 

we  computed  finite  element  solutions  using  uniform  meshes  of  quadratic  triangles. 
In  Table  3.3  (resp.  Table  3.4)  we  give  the  effectivity  indices  for  both  estimators 

in  the  mesh-layers  w,,,  n  =  1,...,  — ,  for  o  =  -  ^resp.  ®  =  g)* 

for  a  =  ^  ^resp.  a  =  the  effectivity  indices  in  the  layers  at  a  fixed  distance 

from  the  singularity  for  both  estimators  converge  to  one  (resp.  zero)  as  the  mesh 
is  refined. 


3.1.S  Uniform  meshes  of  cubic  elements 

We  also  considered  the  mixed  boimdary-value  problem  for  the  Laplace  equar 
tion  in  the  triangular  domain  shown  in  Fig.  3.1c  with  exact  solution 

u(r,  tf)  =  r®  sin(all)  for  a  —  ^  and  a  =  -^  =  ^  ffor  0  =  7  homogeneous 

^  5  \  o 

12 

mixed  boundary  conditions  (resp.  for  a  =  homo^neous  botmdary-conditions) 

o 

were  applied  on  the  edges  Aj  A2,  A|  A3  which  emanate  from  the  singular  point  vhile 
Neumann  botmdary  conditions  consistent  with  the  exact  solution  were  applied  on 

A3A3^  and  we  computed  finite  element  solutions  using  uniform  meshes  of  cubic 

trian^es.  In  Table  3.5  Table  3.^  we  ^ve  the  effectivity  indices  for  the 

mesh-layers  a;„  for  a  =  -^  ^resp.  ®  ~  g)*  W®  observe  that  for  a  =  -g-  ^resp. 

a  =  g J  the  effectivity  indices  in  the  layers  at  a  fixed  distance  from  the  singularity 
converge  to  one  (resp.  zero)  as  the  mesh  is  refined. 

3.1.4  Summary  of  numerical  results 

From  the  results  we  make  the  following  observations.  For  uniform  meshes: 
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'nible  3.S.  Pollution  error  for  uniform  grids.  Convex  domain  shown  in  Fig.  3.1b,  a  =  ~,  p  s  2  (2a  >  p). 

9 

Effectivity  index  for  the  element-residual  and  ZZ  error  estimators  for  each  mesh-layer  Note  that  the 
effectivity  indices  in  the  interior  mesh-layers  are  close  to  one. 


Effectivity  index  for  eadi  mesh-layer 
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Table  S.S.  (continued) 


Effectivity  index  for  each  mesh-layer 

Exact  solution  u(r, 

$)=  rlsinCf 

);  quadratic  elements 

Layer 

h  = 

1 

32 

iiiiiiipcr 

n 

0.647 

0.682 

81.64 

0.995 

1.064 

37.99 

0.997 

0.998 

23.32 

0.976 

0.977 

17.25 

5 

0.967 

0.971 

13.87 

6 

0.964 

0.967 

11.69 

7 

0.962 

0.965 

10.15 

8 

0.960 

0.964 

9.00 

9 

0.959 

0.962 

8.11 

10 

0.957 

7.40 

U 

0.956 

0.959 

6.81 

12 

0.954 

0.957 

6.32 

13 

0.953 

0.955 

5.90 

14 

0.951 

0.954 

5.54 

15 

0.950 

0.952 

5.23 

16 

0.948 

0.950 

4.95 

17 

0.946 

0.948 

4.71 

18 

0.945 

0.947 

4.49 

19 

0.943 

0.944 

4.29 

0.941 
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4.12 

21 

0.939 

0.941 

3.96 

22 

0.937 
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3.81 
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2.90 

32 

0.917 

0.956 

2.83 
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Ikble  8.4.  (continued) 


Effectivity  index  for  endi  mesh-l«yer 

Exact  solution  u(r 

,d)  s  r}tin(^);  quadratic  elements 

Layer 

number 

n 

hs 

T 

32 

lll*klll«_  «/kA 

UPHIl? 

0.764 

0.635 

96.94 

1.186 

0.869 

15.04 

0.548 

7.84 
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5.91 

0.281 

5.05 

6 

0.215 

4.55 

7 

0.170 

4.21 

8 

0.139 

3.96 

9 

0.116 

3.67 

10 

0.099 

3.61 

11 

3.47 

12 

0.075 

3.35 

13 

3.25 

14 

0.059 

3.16 

15 

0.054 

3.08 

16 

0.049 
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3.01 

17 

0.045 

2.94 

18 

0.041 

2.88 

19 

0.037 
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2.83 

20 
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21 

0.032 
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2.73 

22 

2.68 

23 
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24 
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25 
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27 

0.022 

2.50 

28 
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2.48 

29 

0.019 

0.019 

2.46 

0.016 

0.016 

2.44 

31 

0.013 

0.013 

2.43 

32 

0.011 

0.011 

2.42 
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l^ble  S.6.  Pollution  error  for  uniform  grids.  TViangular  domain  shown  in  Fig.  3.1c,  as-,; 

5 

(2a  <  p).  Eifectivity  index  for  the  element-residual  and  ZZ  error  estimaUMrs  for  each  mesh-layer 
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(i) .  The  pollution-effect  depends  on  n  relationship  between  the  exponent  o  of  the 

singularity  and  the  polynomial  degree  p  of  the  elements. 

(ii) .  The  pollution-effect  for  fixed  mesh-sixe  h  and  for  fixed  expmient  a  becomes 

more  significant  as  the  polynomial  degree  p  of  the  elements  is  increased. 

(iii) .  For  elements  of  quadratic  degree  or  higher  the  pollution-effect  in  finite- 

element  solutions  of  the  Lapladan  exists  even  for  convex  domains. 

The  same  conclusions  are  expected  to  hold  for  the  class  of  quasi-imiform  meshes 
which  are  often  used  in  practical  computations. 

3.2  Analysis  of  the  pollution-effect  for  uniform  meshes. 

S.2.1  Preliminaries 

Let  us  assume  that  the  mesh  is  uniform  as  in  the  examples.  (The  analysis  holds 
for  the  broader  class  of  quasi-uniform  meshes.  The  majority  of  the  grids  employed 


in  the  practical  computations  are  quasiuniform), 
solution  of  the  following  boundary-value  problem 

The 

error  = 

u  —  is  the 

“A«a=  12  {12  +»‘t) 

t€T^  c£St 

in 

n 

(3.1a) 

on 

Tn 

(3.16) 

e*=  0 

on 

Tz, 

(3.1c) 

The  variational  form  of  this  problem  was  given  above  in  (2.6)  (the  residual  equa¬ 
tion).  Here  by  J/  we  denote  the  Dirac  function  (line-load)  on  the  edge  e  associated 
with  the  element  r  i.e. 


(3.2) 


where  Jt  was  defined  in  (2.7).  By  the  equilibration  procedure  described  in  Section 
2.2.1  (which  can  be  extended  to  elements  of  any  degree  p)  we  can  rewrite  (the 
right-hand  side  of)  (3.1a)  in  the  form 


=  12  Pr 


(3.3a) 


where 
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<S»r 


(3.») 


where  ere  the  Dirac  fuiicti<m8  analogous  as  before  and  such  that 


(3.4e) 


JprV^O  Wv€V, 


(3.4h) 


as  a  result  ci  the  equilibrations  (here  we  assume  p-d^ree  equilibrations  cf  the 
residuals).  FVmn  (3.4)  it  follows  that  J  pr  —  ^  *nd  hence  there  exists 
such  that 

B,(e„  v)  =  j  f,v  V  r  e  B‘(t)  (3.5) 


Farther  we  have 


I  E  B,(e„«)|  , 


where  ir  denotes  the  exact  solution  of  the  equilibrated  element  residual  problem 
(3.5). 


Let  us  now  define  W-  €  Hb  such  that 

•  p 


-AW;  =  in  n 


dW^  _ 


s  0  on  Fi 


(3.7«) 


(3.75) 


H',  >  0  00  T, 


We  have  by  superposUion 


e*(*)  =  E  W"r(*)  V  <r  €  ft 

r€T^ 


(3.7c) 
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Let  Tq  be  aa  dement  ci  internet  and  deaxit  ^  ^  dementc  in 

m-Uyera  around  Tq,  including  Tq 

*'->(t.)  {t  €  r»  I  TC<‘(r.)},  {  n>  }  (3.9.) 

where 

K‘M:=  U  (  U  Af,),  ;VJ(r.)  =  r.,  :=  (J  (3») 

r'€a<— »(r,)  XeN(T>)  XeN(r^ 

Here  N(r)  is  the  set  of  the  three  vertices  of  triangle  r.  We  let 

and  by  superposition 

V,‘-»(T.)  +  V!<“>(t,)  =  (3.104) 

We  define  the  norms 

:»  IIIK‘">(t.)III,.  .  ♦fV.)  ~  ll|Vi"’(To)lll,.  .  (311.) 

*(To)  =  ll|V^"'(To)  +  »^f’(T.)|||,.  =  IIKIII,.  (3.1U) 

It  is  easy  to  see  that,  firom  the  complementary  variational  principle, 

•;,:=ll|3,.lll,.  >*rK)  (312) 

where  denotes  the  exact  solution  of  the  local  problem  (3.4). 

Remark  S.l.  Here  V/"'*(tq)  is  the  component  of  the  error  which  is  due  to  the  resid- 
uak  in  the  element  Tq  and  m-layers  of  elements  around  it.  Analogously,  is 

the  component  of  the  error  due  to  the  residuab  in  the  elements  outside  the  m-layer 
neighborhood  of  Tq. 
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S.2.2  The  near-field  pollution 

In  general,  the  errw  in  an  element  is  influenced  by  the  reriduals  in  the  entire 
mesh.  However,  when  the  solutkm  is  smo(rth,  and  the  mesh  is  uniform,  the  error 
in  an  eleuMot  Tq  is,  for  all  practical  purposes,  influenced  <Hily  by  the  residuals  in 
To  and  in  a  few  mesh-layers  of  elements  around  Tq.  The  effect  of  the  reriduals  in 
a  few  mesh-layers  of  elements  around  Tq  on  the  error  in  rg  will  be  referred  to  as 
the  near-fieU  pdlulion.  Analogously  the  influence  of  the  reriduak  beyond  a  few 
mesh-layers  of  elements  around  Tq  on  the  error  in  rg  will  be  referred  to  as  the  far- 
field  polhtion.  Below  we  show  results  which  indicate  that  ike  near-fieU  pollution 
is  fully  (^counted  for  by  robust  local  error  estimators. 

We  considered  various  meshes  shown  in  Fig.  3.2a-3.2d,  in  which  the  element  Xg 
is  shaded  black  and  a  few  mesh-layers  around  Xg  are  shaded  pay.  We  computed  the 
norms  defined  in  (3.11a)  first  by  considering  only  the  residuals  in  the  element  Xg  and 
then,  in  progression,  the  reriduals  in  all  the  elements  in  4^”*)(xg )  for  m  =  1, 2, 3, . . .. 
We  computed  the  values  of  the  pollution  factor  for  the  energy: 


lim'”(To)IIP+lllV!/’'(To)|l|» 


(3.13) 


In  Table  3.7a  we  pve  the  values  of  for  p  s  1,  2,  3  for  the  Regular 

mesh-pattem  shown  in  Fig.  3.2a.  In  Table  3.7b  we  pve  the  values  of  for 

p—  1,  2,  3  for  the  Criss-Cross  mesh-pattem  shown  in  Fig.  3.2b.  In  Table  3.7c  we 
give  the  results  for  p  s  1,  2,  3  for  the  mesh-pattem  shown  in  Fig.  3.2c  and  in  Table 
3.7d  we  give  the  values  of  for  p  —  1,  2,  3  for  the  mesh-pattem  shown  in 

Fig.  3.2d.  We  also  see  that  9i°\xg)  is  significantly  smaller  than  the  exact  values 
ilkallU'  factor  converges  to  100%  when  the  reriduals  in 

two  layers  of  elements  around  Xg  are  considered.  Also  note  that  as  the  polynomial 
degree  p  of  the  elements  is  increased,  the  ratio  is  practically  equal  to  100% 

for  m  >  1. 

Hence  we  have  ^i^^Xg)  »  |||ej^|||,^  for  the  mesh-pattems  considered  and  since 
7^  llleAlIU  ^  ^  7^  IllefcllU  for  constants  close  to  1,  (sec  Section 

2.3)  we  have  »  ll|Vi^^^(rb)||lr,?  for  both  the  ER-estimator  with  equilibrated 
data  and  the  ZZ-estimator.  Moreover  since  |||e/^|||r,  =  |||I^^"’^(rg)  +  I^’"^(rb)lllT, 
we  have  «  IllcfcllU  only  if  ll|V^"*^lllr,  <  Hence  expresses 

the  far-field  pollution,  or,  in  other  words  is  an  estimate  of  Hju^  —  (u  —  V2^"‘^)|||t|j 

i.e.  it  estimates  how  close  is  to  u  —  (where  m  =  2  for  the  mesh-pattems 
considered  here). 

From  the  numerical  results  of  this  Section  we  can  infer  that:  When  the  solution 
is  smooth  and  the  mesh  is  periodic  we  have 
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(a) 


Figure  3.2.  Uniform  grids  of  different  puttems  with  the  element  shaded  black  and  a 
few  mesh-layers  around  Tq  shaded  in  different  shades  of  gray,  (a)  Regular 
mesh-pattem;  (b)  Criss-Cross  mesh-pattem;  (c)  Special  mesh-pattem  1; 
(d)  Spedal-mesh-pattem  2. 
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Figure  3.2.  (continued) 
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Tible  S.7a.  Effect  of  the  near  field  pollution:  Pollution  factor  for  the  energy  of  the 
error  in  the  element  Tq  due  to  the  residuals  in  m  mesh-layers  surrounding  Tq  for  the 
mesh  shown  in  Fig.  3.2a.  Note  that  the  percentage  of  the  energy  of  the  ent^  in  Tq 
^proaches  100%  when  the  residuals  in  the  first  layer  surrounding  Tq  are  included. 


Pollution  factor  for  the  energy  of  the  error  in 
the  element  Tq  due  to  the  residuals 
in  m  mesh-layers  around  Tq  . 


Regular  mesh-pattem 


Number  of  Layers 
m 


>^"’(ro)IIP  +  lim'"’(n.)lll’ 


X  100 


p=  1 

■s 

2 

B 

3 

0 

66.92 

74.87 

76.64 

1 

84.44 

97.60 

99.99 

Table  S.7b.  Effect  of  the  near  field  pollution:  Pollution  factor  for  the  enagy  of  the 
error  in  the  element  Tq  due  to  the  residuals  in  m  mesh-layers  surrounding  Tq  for  the 
mesh  shown  in  Fig.  3.2b.  Note  that  the  percentage  of  the  energy  of  the  error  in  Tq 
approaches  100%  when  the  residuals  in  two  layers  surrounding  Tq  are  included. 


'Dible  8.7c.  Effect  of  the  near  field  pollution:  Pollution  factor  for  the  energy  of  the 
error  in  the  element  Tq  due  to  the  residuals  in  m  mesh-layers  surrounding  Tq  for  the 
mesh  shown  in  Fig.  3.2c.  Note  that  the  percentage  of  the  energy  of  the  error  in  Tq 
approaches  100%  when  the  residuals  in  two  layers  surrounding  Tq  are  included. 


Pollution  factor  for  the  energy  of  the  error  in 
the  element  Tq  due  to  the  residuals 
in  m  mesh-layers  around  Tq  . 


Special  mesh-pattem  shown  in  Fig.  3.2c 


Number  of  Layers 


66.73 

92.00 

99.99 

99.99 


TWble  S.7d.  Effect  of  the  near  field  pollution:  Pollution  factor  for  the  energy  of  the 
error  in  the  element  Tq  due  to  the  residuals  in  m  mesh-layers  surrounding  Tq  for  the 
mesh  shown  in  Fig.  3.2d.  Note  that  the  percentage  of  the  energy  of  the  error  in  Tq 
^proaches  100%  when  the  residuals  in  two  layers  surrounding  Tq  are  included. 


Pollution  factor  for  the  energy  of  the  error  in 
the  element  Tq  due  to  the  residuab 
in  m  mesh-layers  around  Tq  . 

Special  mesh-pattem  shown  in  Fig.  S.2d 

Number  of  Layers 

m 

lllW"’(To)lll’  „  JOB 

P=  1 

■a 

■SI 

0 

74.56 

83.04 

1 

93.89 

98.36 

2 

99.98 

99.99 

3 

99.99 
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for  m  >1710 


for  every  interior  element  Tq  €  where  nio  denotes  the  numbers  of  layers  of 
elements  around  Tq  (including  Tq)  for  which  the  element-residuals  are  responsible 
(for  all  practical  purposes)  for  the  error  in  Tq  (for  the  meshes  considered  above 
it  is  sufiBdent  to  take  mg  »  2).  In  other  wor^  when  the  solution  is  smooth  and 
the  mesh  is  periodic,  an  error  estimator  takes  into  account  the  residuals  in  all 
mesh-layers  around  any  interior  element. 

So  far  we  addressed  the  energy-norm  |)|  •  |||^^  in  an  element  rg;  of  course  the 
result  will  hold  if  we  replace  Tg  by  a  patch  of  elements  and  replace  #^”*)(rg )  by 
^(’")(u;*)  with  the  obvious  meaning. 

3.2.3  Asjrmptotics  of  the  pollution-effect  for  uniform  meshes 

Let  us  now  analyze  the  asymptotics  of  the  pollution.  Given  ss  €  H  we  can  write 

-  -  (Tb)(»)  =  13  JprGiix),  »  =  1, 2  (3.14o) 

-^('b)(*)=  E  PrGM,  i  =  l,2  (3.U6) 

where  Gj(a;)  is  the  x, -derivative  of  the  Green’s  function  and  the  integrals  have  to 
be  understood  in  the  proper  sense.  We  have  for  m  >  0 

|Z?"‘G.(*)|<^  on  n(*;i2):=  ft-0(*;i2)  (3.15) 

provided  that  fl  is  a  rectangle  and  either  =  0  or  =  0.  (For  general 
polygonal  domains  (3.15)  holds  for  m  <  m,  where  m  depends  on  the  angles  and 
the  type  of  boundary-conditions.) 

Assuming  that  the  exact  solution  u  is  smooth  we  have 

C.*’<(E-lJ)*SC,*’'  (3.16) 

Assume  now  that  (J  f  D  0(x,  B)  and  using  (3.4)  we  have 

Q|r(m) 

l-^(Tb)(*)|  =  1  E  j  PrGi{x)\  =  \  Y.  /»r(G^.(*)-X)l  (3.17) 

r€a(")(ro)''’'  r€*(")(ro) 
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for  every  X€Str,.  Now  using  (3.15)  we  see  that  we  can  select  Xo  ^ 

Tq  we  have 

I  (G!-Xo)^0  (3.18«) 

and 

|^G<(.)-X,|<C^  (3.18*) 

Hence 

o^m) 

l^(r.)(.)|  =  |  §  -X.)|  = 

*'*»  *■**]*  ^ 


=  1  E  Br(*..Gi(*)-Xo)|<(  E  (3.19) 


rS*(«)(r,) 


and  hence 


^»(f"’('b)||  <c^ 


(3.20) 


Because  Hle^HI  >  we  see  that  if  we  select  R  >  the  pollution  is  n^lipble 
asymptotically. 

In  the  above  argument  we  used  the  asstunption  that  the  solution  u  was  smooth 

land  hence  ( ^  17^)^  < Ch^].  If  the  domain  is  convex  and  either  or  Fp  is 
^  r€T^  ' 

empty  then  (3.16)  holds  for  p  =  1.  Assume  now  that  the  domain  is  not  convex  or 
that  we  have  mixed  boimdary-conditions.  Then  we  have 


81^”’,  ,,  .1  ,  Ch^ 
-^(ro)(»)|  <  . 


i=  1,  2 


(3.21) 


and  hence 


IIIK/"’(’b)llk  < 


(3.22) 


Because 
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IIKIIU,  > 


(3.23) 


ti«  Me  that  if  a  <  ~  then  the  pollutka  may  be  significant. 

2p 

Remark  S.2.  In  the  proof  of  (3.22)  we  employed  the  fact  that  the  Green’s  function 
G(c)  and  its  derivatives  Gj(sB),  i  s  1,  2,  are  ringular  in  the  neighborhood  of  the 
vertices  Ai  (see  Fig.  2.1).  B^um  of  this  singular  behavior  the  exponent  (2a  + 1) 
in  (3.22)  is  in<tependent  of  p.  We  have  proven  in  (3.21)  the  upper  estimate.  It  is 
possible  to  prove  similarly,  as  in  the  themem  13.1  of  [65],  that  the  estimate  (3.21) 
is  sharp  i.e.  that  the  coefficient  2o  cannot  be  improved,  in  general. 

Remark  S.S.  The  estimation  of  in  a  non-asymptotic  way  can  be  made  only 
usic^  on  arpostericm  method.  This  will  be  addressed  in  a  fortbcoming  pi4>er. 


S.2.4  Discussion  of  the  numerical  results  in  the  context  of  the  asymp¬ 
totic  analysis. 

In  Section  3.1  we  gave  numerical  results  which  indicate  the  effect  of  the  pollution 
error  on  the  local  quality  of  the  error  estimators.  We  now  discuss  the  results  of 
Section  3.1  based  on  the  asymptotic  analysis  given  above. 

We  define  the  pollution  factors  (or  a  mesh-cell 


ll|vj“’(.^)IIL.  ...» 

+  IIIKf’V)llli. 


(3.24) 


and 


^  III»(!"*(u>4)IIU 
||H^->(u^)||li. 


100% 


(3.25) 


Obviously  if  pollution  in  u;*  is 

negligible  if  and  are  small  (more  precisely  if  — ►  0  as  fc  — ►  0). 

According  to  the  asymptotic  analysis  the  effect  of  the  pollution  is  negligible 
asymptotically  only  when  2a  >  p.  On  the  other  hand  when  2a  <  p  the  pollution 
error  is  the  significant  part  of  the  total  error  (asymptotically)  in  the  interior  of  the 
mesh;  in  that  case  the  ^Jectivity  index  for  any  estimator  (which  is  based  on  local 
computations),  for  elements  in  the  interior  of  the  mesh,  converges  to  zero. 

For  the  example  problems  given  in  Section  3.1  we  considered  interior  mesh-cells 
and  we  computed  the  components  V^"*^(wo)  ** 

their  energies  |||V^"*^(w§)l||^,  |||Vi"*H<4)|||^.  Let  Go  denote  a  fixed  subdomain 
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which  includes  u;{^  in  its  interkv;  we  mny  nssume  thnt  the  mesh  is  such  thnt  (Ig 
is  always  exactly  covered  by  the  elements  which  intersect  it.  We  let 

Vx  (no)(«)  :=  E  /  Pr  G(m) ,  V,  T,  J  Pr  G(m)  (3.26) 

rMV  r«T. 

'CO,  »«o. 

Below,  for  simplidty,  we  will  omit  Og  fiK»n  most  i^pearanoes  of  ,  V, .  We  note 
that  Vi  and  satisfy  the  foUowing  variational  problems: 

(a)  Find  and 


BaiVr , «)  *  E  (3.27a) 

'CO, 

(b)  Find  V,  €  such  that 

Bo(V,,t;)=  E  FM  (3.276) 

'Sn-Q, 

We  note  that  for  mesh-cells  C  Og  which  are  separated  by  two  (or  more) 
mesh-layers  from  the  boundary  of  fig  we  have 

»  V,  «  V,  (n.)|^  (3.28) 

Thtw,  the  fimctioiu  V\  and  axe  practically  the  same  aa 
for  all  mesh-cells  b;Q  which  are  separated  by  two  or  more  mesh-layers  from  the 
boundary  of  the  subdomain  Bg.  Based  on  the  analysis  given  above  we  say  that 
the  pollution  is  significant  if  the  ratio  is  close  to  100%.  (The  pollution-factor 

is  defined  as  in  (3.24)  where  Vj  (Hg)!^,  Vj  (^)|^  may  be  employed  instead  of 

Vi"*^(u)^),  V2"*^(u)ig).)  We  will  refer  to  pj,  as  the  percentage  of  pollution  in  the 
mesh-cell  Uq. 

We  considered  a  shrinking  eight  element  mesh-cell  u;o  C  Bg  which  consists  of 
the  eight  elements  in  the  A-neighborhood  of  the  center  of  the  subdomain  fig 

5(*°,  A)  :=  |x  =  (x„  Xj)!  max(|x,  -  x?| ,  lx,  -  x5|)  <  a|  (3.29) 

(Note  that  in  all  the  examples  below  we  consider  square  subdomains  fig).  We 
computed  and  in  two  different  ways: 
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а.  Leastsfuttreg  jU  (smoothing) 

We  let  s  u  —  to  and  ,  where  to  is  the  limit  function  to  which 

converges  in  the  neighborhood  of  as  h  tends  to  zero.  The  fimctitm  to  is 
approxiinated  by  letting 

«>(*)  «  o^  +  (a{  ^(*  -  *®)  +  ot}  Qi(m  -  *®))  ,  (3.30o) 

kml 

where 

(?}(«)  =  7le(z*) ,  Q5(*)  =  ,  z  =  *,  +  tij ,  *  =  (xi.xa)  (3.306) 

and  by  employing  a  least-squares  fit  of  the  above  expression  to  the  values  of  U|^  in 
flo  to  determine  the  coefficients. 

б.  Approximation  ofVi,V^  using  a  global  hi^er^order  approximation 

We  can  iq>proximate  by  solving  directly  the  residual  equations  (3.28a) 

(resp.  (3.28b))  satisfied  by  (resp.  )  using  a  finite-element  approximation  of 
degree  g  >  p  on  the  same  mesh  used  to  approximate  U|^  (we  call  this  q-degree 
^proximation  an  overkill). 


In  the  results  below  we  report  the  vdiues  of  the  pollution  factors  the 

effectivity  index  «  ^  and  the  effectivity  index  k\  ;=  —  for  the  shrinking 

^  ^  IIIm 

eight-element  mesh-cell  located  at  the  center  of  the  subdomain  Qq. 

S.2.4’1-  Uniform  meshes  of  linear  elements. 


In  Tables  3.8a,  3.8b  we  give  the  values  of  j||Vj  |||^,  ||m  |||^  and  Wj  •  Wj , 
the  values  of  the  pollution  factors  and  and  the  effectivity  indices  and 

2 

for  the  two  estimators  for  linear  elements  (p  —  1)  and  oc  =  -z  (2a  >  p).  Note 

®  V 

that  the  pollution  factors  and  converge  to  zero  and  the  effectivity  indices 
Kj,  and  Kjt  (for  both  estimators)  converge  to  one  as  the  mesh-size  tends  to  zero. 
2 

Hence  for  a  =  ~  there  is  no  pollution  for  linear  elements,  asymptotically.  In  Fig.  3.3 

(resp.  Fig.  3.4)  we  show  the  directional  derivatives  of  the  V|  -  and  -components 
of  the  error  along  the  s-axis  shown  in  Fig.  3.1a,  for  h  =  0.1^  (resp.  h  =  0.03125). 
In  this  example  (and  all  other  examples  involving  the  ^-shaped  domain)  we  let 
ft  =  (—1,1)*  —  [0,1]  X  [—1,0].  Here  we  let  ftg  =  (0.25,0.75)*.  The  points  of  the 
graph  which  lie  in  the  interior  of  the  subdomain  ftg  are  between  the  dotted  vertical 
lines.  Note  that  the  V,  -component  of  the  error  is  practically  zero  in  the  interior  of 
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Tkble  S.8b.  PoUuti<m  facUMra  fiu^%  and  C.^%  and  effectivity  indices  in  the  shrinking 
mesh-cdl  u^.  L-shi4>ed  domain  shown  in  Fig.  3.1a,  a  s  |  p  s  1  (2a  >  p).  Note  that 
the  pollution  factors  amverges  to  zero  and  the  effectivity  indices  and  amvetge 
to  one,  asymptotically. 


Pollution  factors  and  effectivity  indices  in 


Exact  solution  u(r,8)  s=  sin(^);  linear  elements 


Shrinldng  mesh>cell  centered  at  (0.5,  0.5) 


■  TiMt  inrt 


40.51 

44.31 

30.96 

32.56 

23.19 

23.84 

16.28 

16.50 

Directional  Derivative  of  the  Error 


0.100 


Study  of  the  Pollution  Error 
Singular  Solution:  u*x(r,0)=  Sin(26/3) 

Lineair  Elements.  Uniform  Mesh!  (h  0.125) 


0.050  4 


-0.000 


-0.050  H 


V,  Error  From  Overkill  (p  =  3) 

V,  Error  From  Overkill  (p  =  3) 

Vj  Error  Frjom  Smoothing  (p  =  6) 

Vg  Error  Frpm  Smoothing  (p  =  6) 


■0.100 

0. 


18 


TT 


T  I  I  I  I  ■  I  I 


I  I  I  l“T  I' 


Arc 


0.58  0.98  1.38 

Length  Coordinates  Along  S-Axis 


Figure  3.3.  Pollution  effect  for  uniform  meshes.  Linear  dements  (p  =s  1),  a  =  -, 

1  ^ 
Uniform  mesh  nze  A  =  -.  Directional  derivatives  of  the  -  and  Vg  - 


components  of  the  error  along  the  line  A,  A.  (shown  in  Fig.  3.1a).  Note 

I  OYr  I  •  1  a  «  \  / 


that 


dVj 
*  ^ 

dv. 

ds  ^ 

ds 

in  the  interior  of  the  subdomain. 
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Directional  Derivative  of  the  Error 


0.050 


0.030 


0.010 


Study  of  the  Pollution  Error 

Singular  Solution:  Ua(r.6)  =  Sin(2@/3) 

linear !  Elements,  Uniform  Mesh  (U  =  0.03125) 


-0.010 


-0.030 


-0.050  +-P 
0.14 


••••«  Vj  Error  From  Ovprkill  (p  h  3) 

Va  Error  From  Ovprldll  (p  =  3) 

«-»«-•«  Error  From  Sniioothing  (p  =  6) 

——M  Vj  Error  From  Snioothing  (p  =  6) 

ZZ  Error  Estimatcjr 


0.54  0.94 

Arc  Length  Coordinates  Along  S-Axis 


Figure  3.4.  Pollution  effect  for  uniform  meshes.  Linear  dements  (p  =  1),  a  =  -, 

1  ^ 
Unifonn  mesh  sixe  h  =  — .  (a)  I^rectional  derivatives  of  the  -  and  - 

components  of  the  error  and  the  ZZ  estimators  along  the  s-aads  (shown 

in  Fig.  3.1a).  (The  extent  of  the  subdomain  is  indicated  by  the  dotted 

vertical  lines.)  (b)  Detail  of  the  graph  in  Fig.  3.4a  for  the  pmnts  lying  in 

the  subdomain.  Note  that  and  {tr*^  —  V«il)  •  a 

I  os  I  I  vs  I  os 

asymptotically,  in  the  interior  of  the  subdomain. 
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Directional  Derivative  of  the  Error 


0.010 


Study  of  the  Pollution  Error  ^  ^ 

Singular  Solution:  un(r.6)  =  Sin(2@/3) 

linear  Elements.  Uniform  Mesh  (h  —  0.03125) 


0.005 


0.000 


-0,005  ^ 


-0.010 


0.53 


Vi  Error  From  Overkill  (p  =  3) 
Va  Error  From  Overkill  (p  =  3) 
V,  Error  From  Smoothing  (p  = 
Vj  Error  From  Smoothing  (p  = 
ZZ  Error  Estimator 


I  I"  I  I  ■cr 

0.63 


I  I  I  -I  t  I  I  t  I  I  f  I  t  r  r-r 


■T  i"r~r 


0.73 


0.83 


Arc  Length  Coordinates  Along  S-Axis 


6) 

6) 


(b) 


Figure  3.4.  (continued) 


65 


the  subdomain.  In  Fig.  3.4b  we  show  the  detiul  of  the  graph  in  Fig.  3.4a  for  the 
elements  in  the  interior  of  the  subdomain.  Note  that  the  error  estimator  coincides 
with  the  Vj  -component  of  the  error  in  the  interior  (one  layer  of  elemoits  away 
from  the  boundary)  of  the  subdomain  fig. 

In  Tables  3.9a,  3.9b,  we  give  the  values  of  j||Vi  |||^,  |||  V,  |||^  and  Wj  •  W, , 

the  values  of  the  pollution  factors  ft^  and  and  the  effectivity  indices  and 

for  the  two  estimators  for  linear  elements  (p  =  1)  and  a  ~  r  (^o  <  p).  Note 
"  o 

that  the  pollution  factor  converges  to  100%,  the  pollution  factor  diverges  to 

infinity  and  the  effectivity  index  ic  ^  converges  to  zero;  however  the  effectivity  index 

“0 

k'^  converges  to  one  as  the  mesh-size  tends  to  zero.  In  Fig.  3.5  (resp.  Fig.  3.6)  we 
show  the  directional  derivatives  of  the  components  of  the  error  along  that  part  of 
AjAs  which  lies  over  the  subdomain,  for  h  =  0.125  (resp.  0.03125).  Note  that 
the  V2  -component  of  the  error  converges  to  zero  slower  than  the  -component  as 
the  mesh  is  refined.  In  Fig.  3.6b  we  show  the  detail  of  the  graph  in  Fig.  3.6a  for 
elements  in  the  interior  of  the  subdomiun.  Note  that  the  error  estimator  coincides 
with  the  Vi  -component  of  the  error  in  the  interior  of  the  subdomain. 


S.2.4.2.  Uniform  meshes  of  quadratic  elements. 

Similarly  as  in  the  example  3.1.2  we  let  ft  be  the  trapezoid  shown  in  Fig.  3.1b 
with  Aj  =  (0, 0),  Aj  s:  (1, 0),  A3  =  (1, 1),  A4  =  (-1, 1).  For  the  .«;finition  of  VJ  and 
V2  and  the  shrinking  mesh-cell  we  select^  the  subdomain  ftj.  =  (0.375,0.625)^. 
In  Tables  3.10a,  3.10b  we  report  the  values  of  the  components  of  the  error,  the 
l>ollution  factors  and  the  effectivity  indices  for  quadratic  elements  (p  s  2)  and 

a  =  -  (2a  >  p).  Note  that  the  pollution  factors  and  converge  to  zero 
and  the  effectivity  indices  and  for  both  estimators  converge  to  one  as  the 
mesh-size  tends  to  zero.  In  Tables  3.11a,  3.11b  we  give  the  values  of  |||Vi  |||^, 
lim  lll^  and  fjk  Wj  •  VV^  ,  the  values  of  the  pollution  factors  and  and 

2 

the  effectivity  indices  kk  and  k'k  for  quadratic  elements  and  a  =  ;-.  Note  that 

2  *  *  -  ^ 
for  a  =  -  (2a  <  p),  the  pollution  factor  converges  to  100%,  the  pollution 

factor  A  diverges  to  infinity,  the  effectivity  index  a  »  converges  to  zero;  however 
the  effectivity  index  converges  to  one  as  the  mesh-size  tends  to  zero. 

In  Fig.  3.7a  (resp.  Fig.  3.7b)  we  show  the  directional  derivatives  of  the  - 
and  Vj  -  components  of  the  error  along  that  part  of  the  s-axis  which  lies  over  the 

2 

subdomain  ftg  of  the  uniform  mesh  of  quadratic  elements  for  a  =  -  and  h  —  0.0625 

(resp.  h  =  0.03125)  as  shown  in  Fig.  3.1b.  Note  that  the  -component  of  the  error 
in  the  interior  of  the  subdomun  converges  to  zero  slower  than  the  Vj  -component  of 
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l^ble  S.9b.  Pollution  factors  fi^y  and  effectivity  indices  #c^  in  the  shrinking 

mesh-cell  ujq.  Xfshaped  domain  shown  in  Fig.  3.1a,  a  =  p  =  1  (2a  <  p).  Note  that 

o 

the  pollution  factor  converges  to  100%,  the  pollution  factor  diverges  to  infinity 
and  effectivity  index  ic  ^  converges  to  zero  asymptotically.  However  the  effectivity 
index  converges  to  one. 


Pollution  factors  and  effectivity  indices 

Exact  solution  u(r,d)  s  Bind)*  linear  elements 

Shrinking  mesh-cell  centered  at  (0.5,  0.5) 
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Directional  Derivative  of  the  Error 
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Figure  3.5.  Pollution  effect  for  uniform  meshes.  Linear  dements  (p  =  l)t  a  =  -, 

1 

Uniform  mesh  size  h  =  -.  Directional  derivatives  of  the  -  and  - 

8 

components  of  the  error  and  the  ZZ  estimator  in  the  subdomain  along  the 


s-axis  shown  in  Fig.  3.1a.  Note  that 

dv^ 


ds 


l?51 

d» 


and  •  s 


asjrmptotically,  in  the  interior  of  the  subdomain. 
os 
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Directional  Derivative  of  the  Error 
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Study  of  the  Pollution  Error 

Singular  Solution:  u^r.O)  =  r  Sin(©/3) 
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6) 
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(a) 


Figure  3.6.  Pollution  effect  for  uniform  meshes.  Linear  dements  (p  =  1),  a  =  -, 

1  ^ 
Uniform  mesh-size  h  =  — .  (a)  Directional  derivatives  of  the  \\  -  and  - 

components  and  the  ZZ  estimator  in  the  subdomain  along  the  s-axis  in 

Fig.  3.1a.  (b)  Detail  of  the  graph  in  Fig.  3.6a  for  dements  in  the  interior 


of  the  subdomain.  Note  that 


avj 

av; 

ds  ^ 

as 

and  -  VuJ  •  a  « 


asymptotically,  in  the  interior  of  the  subdomain. 
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Directional  Derivative  of  the  Error 
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Study  of  the  Pollution  Error 

Singular  Solution:  U|x(r.0)  =  r'*^*  Sin(6/3) 
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Figure  3.6.  (continued) 
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S.lOa.  Components  of  the  energy-norm  of  the  error  in  Convex  domain  shown  in  Fig.  3.1b, 
p  —  2  (2a  >  p).  The  Vj  and  V,  components  of  the  error  were  determined  from  the  sdution  for 


2 


Ikble  S.lOb.  Pollution  factors  ^  and  effectivity  indices  isf^  in  the  shrinking 

mesh-cell  u^.  Convex  domain  shown  in  Fig.  3.1b,  a  s  1  p  =  2  (2a  >  p).  Note  that 

o 

the  pollution  factors  decreases  and  the  effectivity  indices  and  converge  to  one 
asymptotically. 


Pollution  factors  and  effectivity  indices 

Uniform  mesh;  Shrinking  mesh-cell  centered  at  (0.5,  0.5) 

Exact  solution  u(r,0)  =  r^8in(y);  quadratic  elements 
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Directional  Derivative  of  the  Error 


0.0005 


Study  of  the  Pollution  Error 

Singular  Solution:  U|^r,6)  =  Sin(26/3) 
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Figure  3.7.  Pollution  effect  for  uniform  meshes.  Quadratic  elements  p  =  2,  a  =  x- 

(a)  Uniform  mesh-size  k  =  — .  (b)  Uniform  mesh-size  h  =  — .  Directional 

derivatives  of  the  -  and  -components  of  the  error  and  the  ZZ  error 

estimator  in  the  subdomain  along  the  s-axis  shown  in  Fig.  3.1b.  Note  that 

^^1  >  “5^1  and  (<r^^  —  •  s  w  “yniplohcally,  in  the  interior 

as  I  os  I  os 

of  the  subdomain. 
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Directional  Derivative  of  the  Error 


Study  of  the  Pollution  Error 
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Figure  3.7.  (continued) 
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Ikble  S.llb.  Pollution  factors  and  effectivity  indices  in  the  shrinking 

.  2 

mesh-cell  Mq.  Ccmvex  domain  shown  in  Fig.  3.1b,  a  =  p  =  2  (2a  <  p).  Note  that  the 

o 

pollution  factor  conv^ges  to  100%,  the  pollution  factor  diverges  to  infinity,  and 
effectivity  index  converges  to  zero  asymptotically.  However  the  effectivity  index 
converges  to  one. 
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the  error  as  the  mesh  is  refined.  Also,  note  that  the  error  estimator  closely  follows 
the  Vj  -otunponent  of  the  error  in  the  interior  of  the  subdomain  Oq' 

S.S.4.S.  Uniform  meshes  of  cubic  dements. 


Similarly  as  the  example  3.1.3  we  let  (I  be  the  triangular  domain  of  Fig.  3.1c 
with  Ai  =s  (0,  0),  A]  =  (0,  1),  A3  =:  (cos  75*,  sin  75*)  and  selected  a  subdomain 
Oo  u  a  parallelogram  with  vertices  at  (0.33,  0.24),  (0.45,  0.24),  (0.48,  0.36),  (0.36, 
0.36).  In  Fig.  3.8a  (resp.  Fig.  3.8b)  we  show  the  directioniil  derivatives  of  the 
V|  -  and  V2  -  components  of  the  error  along  the  s-axis,  as  shown  in  Fig.  3.1c.  We 

considered  a  uniform  mesh  of  cubic  elements  for  o  =  -  (2a  <  p)  and  k  =  0.03125 

o 

(resp.  h  —  0.0078125).  Note  that  the  V3  -component  of  the  error  converges  to  zero 
slower  than  the  -component  as  the  mesh  size  tends  to  zero.  Also  note  that  the 
error  estimator  closely  follows  the  Vj  -component  of  the  error. 


Remark  S.4.  For  uniform  meshes  both  norm  |||Vi  ||K  and  ||m  |||^  converge  to 
zero  as  the  mesh-size  tends  to  zero.  When  2a  >  p  (resp.  2a  <  p)  IIIV3  |||^ 


converges  to  zero  at  a  faster  (resp.  slower)  rate  that  ||m  |||^.  When 


lim  HU 
Ilivi  HU 


0  (resp.  00)  we  have  no  pollution  (resp.  we  have  pollution),  asymptotically. 


Remark  S.5.  The  values  of  the  Vj  -  and  V,  -components  of  the  error  obtained  using 
the  overkill  are  more  accurate  than  the  values  obtained  by  employing  the  least- 
squares  fit  over  the  subdomain.  Nevertheless,  for  mesh-cells  near  the  center  of  the 
subdomain  the  values  of  the  V|  •  and  V,  -functions  computed  by  both  methodologies 
are  practically  the  same. 


3.3  Summary  of  the  results  for  the  pollution-effect  in  uni¬ 
form  meshes 

In  sununary  we  observe  that  the  results  agree  with  the  predictions  of  the 
asymptotic  analysis,  namely; 

(i) .  For  2a  >  p,  the  pollution  factors  in  the  shrinking  mesh-cell  converge  to  zero 

and  the  effectivity  indices  k'^  for  both  estimators  converge  to  one  as 
the  mesh-size  tends  to  zero. 

(ii) .  For  2a  <  p,  the  percentage  of  pollution  in  the  shrinking  mesh-cell  con¬ 

verges  to  100%,  the  pollution-factor  A  diverges  to  infinity  and  the  effectivity 
index  k  .  for  both  estimators  converges  to  zero.  On  the  other  hand  the  ef- 
fectivity  index  converges  to  one. 
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Directional  Derivative  of  the  Errors 
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Study  of  the  Pollution  Error 

Exact  SoluUon:  Un(r.e)  =  r*^*  Sin(6e/5) 

Cubic  Elements;  Uniform  Ifesh  (h  =  0.03125) 


(a) 


Figure  3.8.  Pollution  effect  for  uniform  meshes:  Cubic  element  (p  =  3),  a  = 
Directional  derivatives  of  the  -  and  -components  of  the  error 
the  s-axis  shown  in  Fig.  3.1c.  (a)  Uniform  mesh  size  ^  ~  ^  (^)  Uniform 

mesh  uze  h  =  Note  that  ^  •  a  « 

64  I  ds  I  I  ds  I  '  da 

asymptotically,  in  the  interior  of  the  subdomain. 
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Study  of  the  Pollution  Error.., 
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(b) 


Figure  3.8.  (continued) 


).  The  pollution-effect  in  uniform  meshes  becomes  more  severe  as  the  degree  of 
the  elements  is  increased.  This  can  be  clearly  seen  by  comparing  the  graphs 
in  Figs.  3.8b,  3.9b,  3.10b,  where  the  difference  between  the  derivatives  of 
the  Vj  -  and  Vj  -components  of  the  error  in  the  interior  of  the  subdomain  flo 
increases  with  the  polynomial  degree  of  the  elements. 
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4  Control  of  the  pollution-effect  u*>ing  adaptive 
grids 

In  this  Section  we  will  give  numerical  examples  which  demonstrate  that  when 
the  mesh  is  globally  adaptive  (i.e.  nearly  equilibrated  in  the  energy-norm),  the 
poUution  is  controlled. 


4.1  The  h-adaptive  scheme 

Here  we  summarize  the  A-adaptive  scheme  (see  also  [67],  [68])  which  was 
employed  to  construct  the  adaptive  meshes  in  the  numerical  examples  below. 

Given  a  tolerance  6  for  the  relative  error  in  the  global  energy-norm, 

1.  Choose  an  initial  discretization  of  the  domain  ;  let  T|^  —  T^. 

2.  For  the  mesh  T/^,  compute  the  corresponding  finite  element  solution  U|^  and 
obtain  an  estimate  of  the  energy-norm  of  the  error  in  element  r, 

3.  Check  for  convergence.  If 


S  <l|Klllo 


(4.1) 


stop;  otherwise  proceed  to  the  next  step.  Here  N  denotes  the  number  of 
elements  in  the  mesh  T|^. 

4.  Compute  the  target  error  for  the  optimal  mesh  (using  the  principle  of 
equidistribution  of  error)  namely 


®loTyel  “ 


gllKlllo 

Vn 


(4.2) 


5.  For  each  element  r,  predict  the  optimal  local  mesh-size  from  the  formula 


h 


(4.3) 


Here  is  the  predicted  optimal  mesh-size  for  the  subdomain  within  the  el¬ 
ement  r,  r  is  an  exponent  which  depends  on  the  order  p  of  the  approximation 
and  the  regularity  of  the  solution. 
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6.  Construct  a  smooth  representation  of  the  predicted  optimal  mesh-size  over 
the  mesh  which  now  becomes  the  background  grid  (see  [67])  and  use  it  to 
define  the  mesh-size  function  k  :  f2  — »  R***. 

7.  Use  either  h-refinement  or  h-remeshing  to  construct  a  new  grid  77**  with 
minimal  number  of  elements  AT,  such  that  <  h(x),  x  £  r  €  T^T". 

4.2  Assessment  of  the  pollution-effect  for  adaptive  grids 

We  considered  the  mixed  boundary  value  problem  for  the  Laplace  equation  and 

applied  boundary  conditions  consistent  with  the  exact  solution  u(r,  0)  =  rl  sin(70) 

over  the  square-stii  domain  fl  =  (—1,1)’  —  [0,1]  x  {0}  shown  with  thick  black 

perigram  in  Fig.  4.1a.  (Note  that  in  this  case  0  =  7  and  the  singularity  is  very 

4 

strong.)  We  constructed  nearly  equilibrated  meshes  of  linear,  quadratic  and  cubic 
elements  using  adaptive  refinement  and  remeshing  and  we  studied  the  effect  of  the 
pollution  for  interior  mesh-cells. 

4.2.1  Nearly  equilibrated  grids  of  linear  elements 

In  Fig.  4.1a  and  Fig.  4.1c  we  show  the  adaptive  grids  of  linear  elements  (p  =  1) 
generated  using  refinement  and  remeshing,  respectively  (a  tolerance  of  0.5%  for 
the  global  relative  error  in  the  energy-norm  was  employed).  In  Fig.  4.2  (resp. 
Fig.  4.3)  we  show  the  graph  of  the  values  of  the  energy-norms  of  the  error  in  the 
elements  for  the  mesh  shown  in  Fig.  4.1a  (resp.  4.1c)  (The  elements  have  been 
munbered  according  to  their  distance  from  the  singular  point  and  two  rings  of 
elements  around  the  singularity  have  been  excluded  from  the  graph).  The  values 
of  the  effectivity  index  for  the  mesh-cells  (which  are  shown  shaded  gray  in  the 
Figures)  are  reported  in  Table  4.1  for  the  mesh  shown  in  Fig.  4.1a  and  Table  4.2 
shown  in  Fig.  4.1c.  These  values  are  'Compared  with  the  pollution-free  values  of  the 
’  effectivity  index,  «,  for  the  mesh-cells.  The  pollution-free  value  of  the  effectivity 
index  in  each  mesh-cell  was  obtained  from  a  pollution-free  finite  element  solution 
of  the  boundary  value  problem  in  the  domsun  U,  using  the  same  mesh,  with  data 
obtained  from  the  (p  + 1)  degree  Taylor  series  expansion  of  the  exact  solution  about 
the  center  X  of  each  mesh-cell.  Note  that  the  computed  values  of  the  effectivity 
index  are  very  close  to  the  corresponding  pollution-free  values  in  all  the  interior 
mesh-cells. 

For  the  mesh  shown  in  Fig.  4.1c  we  also  computed  the  effectivity  index  for  the 
mesh-layers,  which  are  shown  in  Figs.  4.4a,  4.4b.  We  excluded  the  last  mesh-layer 
of  elements  adjacent  to  the  singular  point  and  numbered  the  layers  inward  from 
the  boundary  as  shown  in  Fig.  4.4b.  In  Table  4.3  we  give  the  computed  effectivity 
indices  and  pollution-free  value  of  the  effectivity  index  for  each  mesh-layer.  The 
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(a) 

Figure  4.1.  (a)  A  iinite-elemeat  grid  ofUnear  elements  constructed  over  the  slit-domain 
using  adaptive  refinement  and  a  tolerance  of  0.5%  for  the  global  relative 
error  in  the  energy-norm;  (b)  Detail  of  the  finite-element  mesh  shown  in 
Fig.  4.1a  near  the  singular  point  with  the  mesh-cell  shown  shaded 
gray;  (c)  A  finite-element  grid  of  linear  elements  constructed  over  the  slit- 
domain  wing  adaptive  remeshing  and  a  tolerance  of  0.5%  for  the  ^obal 
relative  error  in  the  energy-norm. 
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(b) 


Figure  4.1.  (continued) 
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study  of  the  Pollution  Error 
Singular  Solution:  Ujg^r.Q)  =  r*'  '  Sin(0/4) 
Linear  Elements;  Adaptively  Refined  Mesh 
Equidistribution  of  the  Error  in  the  Elements 


Target  Error 


Element  Number 


Figure  4.2.  Control  of  the  pollution-effect  using  adaptive  grids.  Distribution  of  the 
energy-norm  of  the  error  in  the  elements  of  the  mesh  shown  in  Fig.  4.1a. 
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Energy-Norm  of  Error  in  the  Element 


Element  Number 


Figure  4.3.  Control  of  the  poUution-effect  using  adaptive  grids.  Distribution  of  the 
energy-norm  of  the  error  in  the  elements  of  the  mesh  shown  in  Fig.  4.1c. 
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Tkble  4.1.  Control  of  the  pollution-effect  using  adaptive  grids.  Ad^tive  grid  of 
linear  elements  generated  over  the  slit-domain  using  refinement  (tolerance  ss  0.5%, 
error  control).  Comparison  of  the  effectivity  index,  Kj^  ,  in  the  i7*-n<mn,  computed 
with  respect  to  the  exact  solution  with  the  corresponding  values  of  the  pollution-free 
effectivity  index  for  the  ZZ  estimator  and  the  element-residual  (ER)  estimator  in 
the  mesh-ceUs  u)x^ ,  and  (which  consist  of  the  elements  connected  to  the  nodes 
Xx,  AT}  and  shown  in  Fig.  4.1a). 


Adaptive  grid  using  refinement,  H'*-error  control,  linear  elements 


Effectivity  index  in  the  ^*>norm 
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Ikble  4.2.  Control  of  the  pollution-effect  using  adaptive  grids.  Adaptive  grid  of 
linear  elements  generated  over  the  sUt-domain  using  remeshing  (tolerance  =  0.5%, 
J7*-error  norm  equidistribution).  Comparison  of  the  effertivity  index  Kj^  in  the 
norm,  computed  with  respect  to  the  exact  solution  with  the  corresponding  value  of  the 
pollution-free  effecti vity  index  ,  for  the  ZZ  estimator  and  the  element-residual  (ER) 
estimator  in  the  mesh-cells  end  tujf,  (which  consist  of  the  elements  connected 

to  the  nodes  Xi,  and  Xz  shown  in  Fig.  4.1c). 


Adaptive  grid  using  remeshing,  fl’^-error  control,  linear  elements 

Effectivity  index  in  the  if*-norm 

Mesh-cell 

Mesh-cell 

Mesh-cell 

^zz 

=ZZ 

1^ 

jiER 

< 

M 

=zz 

1^ 

M 

m 

1 

0.95 

0.93 

B 

0.94 

0.93 

B 

g 

g 

89 


(a) 

Figure  4.4.  (a)  The  adaptive  finite  element  mesh  in  Fig.  4.1a  in  which  we  indicate 
the  mesh-layers  for  which  we  computed  the  effectivity  index;  (b)  Detail  of 
the  finite-element  mesh  shoadng  the  mesh  layer  near  the  singular  point;  (c) 
Detail  of  the  mesh  near  the  singular  point.  This  is  the  same  mesh,  shown 
in  Fig.  4.4a  and  4.4b,  with  two  rings  of  dements  around  the  singular  point 
removed.  This  mesh  was  used  in  the  computation  of  a  pollution-free  finite 
dement  solution. 
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(b) 


(c) 


Figure  4.4  (continued) 


91 


Table  4.S.  CSontrol  of  the  pollution-effect  using  adaptive  grids.  Adaptive  grid  of  linear 
elements  shown  in  Fig.  4.4.  Effectivity  index  for  element-residual  and  ZZ  estimators 
for  each  mesh-layer  in  the  interior  of  the  mesh.  Note  that  the  effectivity  index 
computed  with  respect  to  the  exact  solution  is  very  close  to  the  pollution-free  values 
for  all  the  layers  in  the  interior  of  the  mesh. 


Adaptive  grid  using  refinement, 

error  control,  linear  elements 

EflTectivity  index  in  each  mesh  layer 
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pollution-free  values  for  the  mesh-layers  were  obtained  by  solving  the  boundary 
value  problem  over  the  domain  which  excludes  two  rings  of  elements  around 
the  singular  point  as  shown  in  Fig.  4.4c.  Note  that  the  values  of  the  effectivity 
index  for  the  mesh-layers  are  very  close  to  the  corresponding  pollution-free 
values  of  the  effectivity  index. 

4.2.2  Nearly  equilibrated  grids  of  quadratic  elements 

In  Fig.  4.5  we  show  a  grid  of  quadratic  elements  (p  =  2)  which  was  generated 
using  adaptive  refinement  and  equidistribution  of  the  energy-norm  of  the  error 
for  a  tolerance  of  0.1%.  In  Fig.  4.6  we  show  the  energy-norm  of  the  error  in 
the  elements  for  the  mesh  shown  in  Fig.  4.5.  (The  elements  have  been  numbered 
according  to  their  distance  from  the  singular  point  and  two  rings  of  elements  around 
the  singularity  have  been  excluded  from  the  graph.)  In  Table  4.4  (resp.  Table  4.5) 
we  compare  the  effectivity  index  of  the  estimators  in  the  shaded  mesh-cells  (resp. 
mesh-layers)  with  the  pollution-free  values  of  the  effectivity  index  in  the  mesh-cells 
(resp.  mesh-layers).  Note  that  the  values  of  the  effectivity  index  for  the  mesh-cells 
and  the  mesh-layers  aie  very  close  to  the  corresponding  pollution-free  values  of  the 
effectivity  index. 

4.2.3  Nearly  equilibrated  grids  of  cubic  elements 

In  Fig.  4.7a  we  give  a  grid  of  cubic  elements  (p  =  3)  which  was  constructed 
using  adaptive  refinement  and  equidistribution  of  the  energy-norm  of  the  error  for  a 
tolerance  of  0.05%.  In  Fig.  4.8  we  show  the  energy-norm  of  the  error  in  the  elements 
for  the  mesh  shown  in  Fig.  4.7.  (The  elements  have  been  numbered  according  to 
their  distance  from  the  singularity  and  the  first  two  rings  of  elements  around  the 
singular  point  have  been  excluded  from  the  graph.)  In  Table  4.6  (resp.  Table  4.7) 
we  compare  the  effectivity  index  of  the  estimators  in  the  shaded  mesh-cells  (resp. 
mesh-layers)  with  the  pollution-free  values.  Note  that  the  values  of  the  effectivity 
index  for  the  mesh-cells  and  the  mesh-layers  are  very  close  to  the  corresponding 
pollution-free  values  of  the  effectivity  index. 

4.3  Summary  of  the  results  for  the  pollution-effect  in 
adaptive  grids 

In  summary,  we  observe  that  when  the  mesh  is  nearly  equilibrated  in  the 
energy-norm,  the  pollution  in  the  energy-norm  is  negligible  and  the  values  of  the 
effectivity  index  are,  for  all  practical  purposes,  the  same  as  the  pollution-free  values 
for  any  estimator. 
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Figure  4.5.  A  finite-demcnt  grid  of  quadratic  elements  constructed  over  the  slit-domain 
adaptive  refinement  and  a  tolerance  of  0.1%  for  the  global  relative 
error  in  the  energy-norm. 


Energy-Norm  of  Error  in  the  Element 


Figure  4.6.  Control  of  the  pollution-effect  using  adaptive  grids.  Distribution  of  the 
energy-norm  of  the  error  in  the  elements  of  the  mesh  shown  in  Fig.  4.5. 
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Ikble  4.4.  Control  of  the  pollution-effect  using  adaptive  grids.  Ad^tive  grid  of 
quadratic  elements  generated  over  the  slit-domain  using  refinement  (tolerance  =  0.1%, 
-error  norm  equidistribution).  Comparison  of  the  effect!  vity  index  ic^,  computed 
with  respect  to  the  exact  solution  with  the  the  corresponding  values  of  the  pollution- 
free  effect! vity  index  ic^,  for  the  ZZ  estimator  and  the  element-residual  (ER)  estimator, 
in  the  mesh-cells  and  uJx^  (which  consist  of  the  elements  connected  to  the 

nodes  X2  and  X3  shown  in  Fig.  4.5). 
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Tkble  4.5.  Control  of  the  pollution-effect  using  adaptive  grids.  Adi^tive  grid  of 
quadratic  elements  shown  in  Fig.  4.5.  Effectivity  index  for  element-residual  and  ZZ 
estimators  for  each  mesh-layer  in  the  interior  of  the  mesh.  Note  that  the  effectivity 
index  computed  with  respect  to  the  exact  solution  is  very  close  to  the  pollution-free 
values  for  all  the  layers  in  the  interior  of  the  mesh. 


Figure  4.7.  (a)  A  finite-element  grid  of  cubic  elements  constructed  over  the  slit-domain 
iming  adaptive  refinement  and  a  tolerance  of  0.04%  for  the  global  relative 
error  in  the  energy-norm;  (b)  Detail  of  the  finite-element  mesh  shown  in 
Fig.  4.7a  near  the  singular  point  with  the  mesh-cell 
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Figure  4.8.  Control  of  the  pollution-effect  using  adaptive  grids.  Distribution  of  the 
energy-norm  of  the  error  in  the  elements  of  the  mesh  shown  in  Fig.  4.7. 


Table  4.6.  Ck>iitrol  of  the  pollution-effect  using  adaptive  grids.  Adaptive  grid  of 
cubic  elements  generated  over  the  slit-domiun  using  refinement  (tolerance  =  0.05%, 
-error  norm  equidistribution).  Comparison  of  the  effectivity  index  computed 
with  respect  to  the  exact  solution,  with  the  corresponding  values  of  the  pollution-free 
effectivity  index  for  the  ZZ  estimator  the  element-residual  (ER)  estimator,  in  the 
mesh-cells  ci;jf ^ and  (which  consist  of  the  elements  connected  to  the  nodes 
Xu  X%  and  Xz  shown  in  Fig.  4.7a). 
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Tftble  4.7.  Control  of  the  pollution-effect  using  adi4>tive  grids.  Adaptive  grid  of  cubic 
elements  shown  in  Fig.  4.7a.  E^ectivity  index  for  element-residual  and  ZZ  estimators 
for  each  mesh-layer  in  the  interior  of  the  mesh.  Note  that  the  effectivity  index 
computed  with  respect  to  the  exact  solution  is  very  close  to  the  pollution-firee  values 
for  all  the  layers  in  the  interior  of  the  mesh. 


5  Pollution-effect  for  meshes  refined  locally  in 
the  interior 

As  discussed  in  the  Introduction  it  is  a  common  practice  in  engineering  compu¬ 
tations  to  refine  the  mesh  locally,  only  in  local  re^ons  of  interest.  In  this  Section 
we  give  numerical  evidence  which  shows  that:  Local  refinements  tn  a  snidomain 
of  interest  in  the  interior  of  the  mesh  (where  the  solution  is  smooth)  do  not  lead 
to  hi^er  heal  accuracy  unless  the  mesh  is  also  properly  refined  in  Vie  rest  of  the 
domain.  In  particular  we  demonstrate  that  if  the  singular  points  are  not  refined 
and  the  mesh  is  refined  locally  in  an  interior  subdomain 

a.  The  finite  element  solution  in  the  interior  of  the  subdomain  converges  to 
a  fimction  which  satisfies  the  differential  equation  (i.e.  it  is  harmonic)  but 
may  be  very  different  from  the  exact  solution  of  the  boundary  value  problem 
(i.e.  the  finite  element  solution  converges  to  a  wrong  solution  in  the  interior 
of  the  subdomain). 

b.  The  Vi  -component  of  the  error  converges  to  zero  in  the  interior  of  the  sub- 
domain  while  the  -component  remains  practically  constant.  The  error 
estimators  estimate  only  the  \\  -component  of  the  error  in  the  interior  of  the 
subdomain  and  if  ||m  |||q^  >  |||Vj  |||()^  severe  underestimation  of  the  error 
may  occur  in  the  interior  of  the  subdomain. 


5.1  Convergence  of  the  approximate  solution  in  the  inte¬ 
rior  of  the  subdomain 


We  first  present  a  simple  numerical  example  to  demonstrate  that  when  the 
mesh  is  refined  locally  in  an  interior  subdomain  only,  as  shown  for  example  in 
Fig.  5.1,  the  finite  element  solution  converges  to  a  wrong  solution  in  the  interior  of 
the  subdomain. 


We  considered  the  mixed  boundary-value  problem  for  the  Lapladan  in  the 
Xrshaped  domain  with  exact  solution  v(r,d)  =  ri  sin  The  boundary- value 
problem  was  solved  using  the  mesh  shown  in  Fi^.  5.1  where  the  domain  was  meshed 
by  a  coarse  uniform  grid  ^with  mesh-size  hg  s  and  the  subdomain  is  subdivided 

uniformly  five  times  ^the  mesh-size  in  the  subdomain  h  =  ^  hg) ,  and  elements 

of  degree  p  (p  —  1,  2  and  3).  From  the  finite  element  solution  we  approximated 
the  function  w  :=u  —  Vj  in  the  subdomain  and  the  exact  solution  u  by  their  local 
Taylor-series  expansions  about  the  center  of  the  subdomain  x  =  (zj ,  z"),  namely 
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where  t  —  1,2  denotes  the  harmonic  m<m<»nial8  of  d^ree  q;  a _ is  the 

maximum  degree  employed  in  the  harmmuc  Taylor-series  expansions  (bekm  we  let 

a _ =  4);  SE°  =  (0.375,  0.375)  is  the  center  of  the  subdomain. 

In  Table  5.1  we  report  the  ^-coefficients  of  the  local  Taylor-series  expansion  of 
the  exact  solution  u  and  the  corresponding  ot-coeffidents  of  the  local  Taylor-series 
expansion  of  the  limit  function  u;  for  p  =  1,  2  and  3.  Note  that  the  coeffidents  of 
the  limit-function  tv  are  very  different  from  the  coeffidents  of  the  exact  solution  u. 
Hence  if  the  mesh  is  refined  only  in  an  interior  subdomain  and  the  mesh  outside  the 
subdomain  is  kept  fixed  the  finite  element  solution  converges  to  a  wrong  solution 
in  the  interior  of  the  subdomain. 


5.2  Asymptotics  of  the  error  in  the  interior  of  a  localfy 
refined  subdomain 

Here  we  study  the  asymptotics  of  the  components  ((Iq)  and  (Qq)  of 
error  in  the  interior  of  a  subdomain  Hq  as  the  mesh  in  the  subdomain  is  successively 
refined  while  the  mesh  outdde  the  subdomain  remdns  fixed. 


5.2.1  Locally  refined  meshes  of  linear  elements 
5.2.1. 1.  Smooth  solutions 


We  solved  the  boundary-value  problem  for  the  Laplace  equation  in  the  domain 
n  =  (0, 1)^  and  employed  boundary-conditions  consistent  with  the  smooth  har¬ 
monic  solution  u(zj,Z3)  =  itsinxzisinhxZ].  In  Figure  5.2  we  show  the  domain 

12  meshed  with  a  coarse  uniform  grid  with  mesh-siase  Hq  =  -  and  the  subdomain 

1 

12q  :=  (0.375,0.625)^  meshed  with  a  finer  grid  with  mesh-siase  h  sz  We  con- 

o4 

sidered  linear  elements  (p  =  1)  and  different  mesh-sizes  in  the  subdomain  shown 
in  Fig.  5.2.  In  Tables  5.2a,  5.2b  we  give  the  values  of  the  components  |||V|  |||^, 

111^11115*  energy-norm  of  the  error,  the  values  of  the  pol- 

lution  factors  effectivity  indices  and  for  the  shrinking 

eight-element  mesh-cell  at  the  center  of  the  subdomain.  Note  that,  as  the  mesh- 
size  in  the  subdomain  tends  to  zero,  the  pollution  factor  converges  to  100%, 
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l^ible  5.1.  Pollution-effect  for  meshes  refined  locally  in  an  interior  subdomain.  CoefiB- 
dents  of  the  local  Taylor  series  expansions  of  the  exact  solution  u  and  the  limit  function 
w  about  the  center  of  the  subdomain  Qq  (0.5,0.75)^.  Zrshaped  domain  shown  in 
Fig.  5.1,  exact  solution  u(r,0)  =  sin  mesh-size  in  the  subdomain  h  =  Note 
that  the  /9-coeffidents  of  the  local  Taylor  series  expansion  of  u  are  very  different  from 
the  ot-coeffidents  of  the  local  Taylor  series  expansion  of  the  limit  fimction  w. 
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the  pollution  factor  diverges  to  infinity,  the  effectivity  index  converges  to 
aero  and  the  effectivity-index  converges  to  one. 

In  Fig.  5.3a  (resp.  Fig.  5.3b)  we  give  the  graph  of  the  directional  derivative 
of  the  -  and  -components  of  the  error  in  the  subdomain  shown  in  Fig.  5.2 

when  the  mesh-size  in  the  subdomain  is  h  =  (resp.  h  =  )* 

directional  derivative  of  the  V,  -component  of  the  error  is  very  large  relative  to  the 
directional  derivative  of  the  Vj  -component  of  the  error.  As  the  mesh-size  in  the 
subdomain  is  reduced  the  -component  of  the  error  converges  to  zero  while  the 
-component  remains  practically  constant. 

5.2.1. 2.  Singular  solutions 


Next  we  solved  the  mixed  boundary  problem  for  the  Laplace  equation  and 

considered  the  singular  solution  u(r,0)  =&  rt  sin(-tf)  for  linear  elements  (p  s  1) 

and  different  mesh-sizes  in  the  subdommn  shown  in  Fig.  5.1.  In  Tables  5.3a, 

5.3b  we  give  the  values  of  ||m  |||^,  jjlVj  |||^  and  f  •  VV^ ,  the  values  of 

0  0  j 

the  pollution  factors  and  and  the  effectivity  indices  and  for  the 
shrinking  eight-element  mesh-cel!  at  the  center  of  the  subdomain,  as  the  mesh-size 
in  the  subdomain  tends  to  zero.  Note  that,  as  the  mesh-size  in  the  subdomain 
tends  to  zero,  the  pollution  factor  converges  to  100%,  the  pollution  factor 
diverges  to  infinity  for  a  larger  mesh-size  in  the  subdomain  than  in  the  case  of  the 
smooth  solution,  the  effectivity  index  for  both  estimators  converges  to  zero 
and  the  effectivity  index  k\  converges  to  one. 

^  1 

We  also  considered  the  exact  solution  «(r,^)  =  ri  8in(-i9).  In  Fig.  5.4  we  give 
the  graphs  of  the  directional  derivative  of  the  -  and  -components  of  the  error 
in  the  subdomain  shown  in  Fig.  5.1,  along  the  s-axis  shown  in  Fig.  3.1a,  obtained 
by  .employing  an  overkill  with  p  =  3,  when  the  mesh-size  in  the  subdomain  is 

h  =  —hQ.  In  Fig.  5.5  we  show  the  detail  in  the  interior  of  the  subdomain  for  the 
graphs  in  Fig.  5.4.  Note  that  the  directional  derivative  of  the  V2  -component  of 
the  error  is  very  large  relative  to  the  directional  derivative  of  the  -component 
of  the  error  in  the  subdomain.  Also  note  that  the  error  estimator  coincides  with 
the  Vj  -component  of  the  error  in  the  interior  of  the  subdomain.  In  Fig.  5.6  we 
plot  the  directional  derivative  of  the  Vj  -  and  Vj  -components  of  the  error  along 
the  s-axis  shown  in  Fig.  3.1a  when  the  mesh-size  in  the  subdomain  is  reduced  to 

h  =  —  hg.  As  the  mesh-size  in  the  subdomain  is  reduced  the  -component  of 
the  error  converges  to  zero  in  the  subdomain  while  the  -component  of  the  error 
remains  practically  constant. 
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Figure  5.2.  Pollution-effect  for  meshes  refined  locally  in  an  interior  subdomain.  The 

iinit-gquare  domain  meshed  with  a  coarse  umform  gnd  with  mesh  size  Hq  = 

-  The  mesh  has  been  divided  uniformly  five  times  in  an  interior  subdomain 

1  \ 

Hq  (^the  mesh-size  in  the  subdomain  is  h=  2^)* 
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Tftble  5.2a.  Components  of  the  energy-norm  of  the  error  in  Uq.  Square  domain  shown  in  Fig.  5.2,  smooth 
solution  u(x|,X3)  =  AsinirxjsinhirX],  A  =  10^,  p=  1.  The  V^-  and  -components  were  determine  from 
the  solution  for  h  =  using  a  least  square  fit  over  fig  ss  (0.375,  0.625)  x  (0.375,  0.625). 
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l^ble  5.2b.  Pollution-effect  for  meshes  refined  locally  in  an  interior  subdomain. 
Pollution  factors  and  ^  and  effectivity  indices  in  the  shrinking  mesh-cell  u^. 
Square  domain  shown  in  Fig.  5.2,  smooth  solution  u(xi,X2)  =  >4 sin XX|  sinh xx^,  A  — 
1(P,  p  =  1.  Note  that  the  pollution  factor  ocnverges  to  100%,  the  pollution  factor 
diverges  to  infinity  and  the  effectivity  index  converges  to  zero  asymptotically; 
however  the  effectivity  index  k\  conv^ges  to  one  as  the  mesh-size  tends  to  zero. 


Pollution  factors  and  effectivity  indices  in  Ug 


Exact  solution  u(X|,x^)  =  Asinxxjsinhxx,;  linear  elements 
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Directional  Derivative  of  the  Error 


Study  of  the  Pollution  Error 
2.0E-4  -1  Smooth  Solution:  Ub^Xi.Xc)  =  A  SinnxiSinhnxj 

I  Linear  Elements.  Uhiform  Refinements  in  Subdomain  Only 
'Mesh  Size  in  Subdomain  =  0.00390625 
I  Mesh  Size  in  Rest  of  the  Domain  s  0.125 
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Figure  5.3.  (continued) 
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Ikble  S.Sb.  Pollution-effect  for  meshes  refined  locally  in  an  interior  subdomaiu. 
Pollution-factors  and  and  effectivity  indices  in  shrinking  mesh-cell  L- 
s1m4>H  domain  shown  in  Fig.  5.1a,  a  s  p  —  1.  Note  that  the  pollution  factor 
converges  to  100%,  the  pollution  factor  ^  diverges  to  infinity,  and  the  effectivity 
index  converges  to  zero.  However,  the  effectivity  index  converges  to  one. 


Pollution  factors  and  effectivity  indices 


Shrinking  mesh-cell  centered  at  (O.S,  0.5) 
Exact  solution  u(r,<f)  =  rl  8in(^);  linear  elements 


Uniform  mesh-refinement  in  the  subdomain  only 


Mesh-size  in  the  rest  of  the  domain  k  =  0.125 
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Directional  Derivative  of  the  Error 


Figure  5.4.  Pollution-effect  for  meshei  refined  locally  in  an  interior  aubdomain.  Linear 
elements  (p  =  1),  a  =  X-shaped  domain  shown  in  Fig.  5.1.  Uniform 

mesh  size  in  the  subdomain  h  =  Directional  derivatives  of  the  Vj  - 

IJSo 

and  V] -components  of  the  error  obtained  from  the  overktU.  Note  that 

>  1^!]^  and  ^  asymptotically,  in  the  interior 

ds  I  ds  0* 

of  the  subdomain. 
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Direcliorial  Derivative  of  the  Error 
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Study  of  the  Pollution  Error 
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Figure  5.5.  Pollution-effect  for  meshes  refined  locally  in  an  interior  subdomain.  Linear 

elements  (p  =  1),  a  =  L-shaped  domain  shown  in  Fig.  5.1.  Uniform 

^  1 

mesh  size  in  the  subdomain  k  =  Directional  derivatives  of  the  -  and 
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Vg  -components  of  the  error  and  the  ZZ  error  estimator  in  the  subdomain 
Qg.  Comparison  between  the  values  obtained  from  smooihsng  and  overkitt. 
Note  that  the  directional  derivatives  of  andV,  obtained  from  smoothing 

dV 

and  overkill  coincide  exactly  and  ^  ^  ^  in  the  interior  of 

the  subdomain. 


Directional  Derivative  of  the  Error 
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Study  of  the  Pollution  Error 

Singular  Solution:  u^r.3)  =  Sin(6/3) 

Linear  Elements.  Uniform  Refinements  in  Subdomain  Only 
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Figure  5.6.  Pollution  effect  for  grids  refined  locally  in  an  interior  subdomain.  Linear 
elements  (p  =  1),  a  =  Zr8hi4>ed  domain  shown  in  Fig.  5.1,  Uniform 
mesh  size  in  the  subdomain  h  =  Directional  derivatives  of  the  V,  - 

IO¥r  I  /Jl/ 

asymptotically, 

in  the  interior  of  the  subdomain. 
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5.2.2  Locally  refined  meshes  of  quadratic  elements 

5.2.2.1.  Smooth  solution 

We  considered  the  same  boundary-value  problem  as  in  Secti<m  5.2.1.1  with  the 
smooth  harmonic  solution  u(xi,X3)  s  Asinxxisinhxxj  and  employed  meshes  of 
quadratic  elements  (p  —  2)  with  various  mesh-sises  in  the  subdomain  shown  in 
Fig.  5.2.  In  Tbble  5.4a,  5.4b  we  report  the  values  of  the  p<dlution  factors  and  the 
effectivity  indices  for  different  mesh-sises  in  the  subdomain.  Once  again  we  note 
that  the  pollution-factor  converges  to  100%,  the  pollution  factor  ^  diverga 
to  infinity  and  the  effectivity  index  converges  to  xero  as  the  mesh-size  in  the 
subdomain  tends  to  zero. 

5.2.2.2.  Singular  solutions 

We  solved  the  mixed  boundary  value  problem  for  the  Lapladan  with  boundary- 
conditions  consistent  with  the  exact  solution  u(r,0)  s  nn(-tf)  fand  employed 

meshes  of  quadratic  elements  (p  =  2)  with  various  mesh-sizes  in  the  subdomain 
shown  in  Fig.  5.7.  In  Tables  5.5a,  5.5b  we  give  the  values  of  |||Vi  |||^,  ||m  }||^ 

/O  0 

•  VV^ ,  the  values  of  the  pollution  factors  and  and  effectivity 

indices  ic^  and  in  the  shrinking  eight-element  mesh-cell  at  the  center  of  the 
subdomain,  as  the  mesh-size  in  the  subdomain  tends  to  zero.  Again,  we  note 
that  the  percentage  of  pollution  in  the  shrinking  mesh-cell  at  the  center  of  the 
subdomain  converges  to  100%,  the  pollution  factor  diverges  to  infinity,  the 
effectivity  index  for  both  estimators  converges  to  zero  and  the  effectivity  index 
converges  to  one.  Also  note  that  compared  with  the  results  in  Tbble  5.2,  we 
see  that  the  pollution-effect  is  stronger  for  larger  values  of  the  mesh-size  in  the 
subdomain. 

We  also  considered  the  exact  solution  tt(r,9)  =  ri  sin(-tf).  In  Fig.  5.8a  (resp. 

Fig.  5.8b)  we  show  the  graph  of  the  directional  derivative  of  the  -  and  Vj  - 
components  of  the  error  in  the  subdomain  shown  in  Fig.  5.7,  along  the  s-axis  shown 

in  Fig.  3.1b,  where  the  mesh-size  in  the  subdomain  ^  ~  ^  (resp.  h  = 

Note  that  the  directional  derivative  V2  -component  is  significantly  greater  than  me 
directional  derivative  -component  of  the  error  in  the  interior  of  the  subdomidn. 
F\irther  note  that  as  the  mesh-size  in  the  subdomain  is  reduced,  the  Vj  -component 
of  the  error  converges  to  zero  while  the  -component  remains  constant.  Also  note 
that  the  error  estimator  coincides  with  the  -component  of  the  error. 


118 


Table  5.4a.  Components  of  the  energy-norm  of  the  error  in  Square  domain  shown  in  Fig.  5.2,  smooth 
solution  u(X|, Xa)  -  i4sin irz^  sinh  wxj,  A  =  10^,  p-2.  The  -  and  V,  -components  were  determine  frmn 
the  solution  for  h  —  using  a  least  square  fit  over  Hq  s  (0.375,  0.625)  x  (0.375,  0.625). 


119 


Ikble  5.4b.  Pollution-effect  for  meshes  refined  locally  in  an  interior  subdomain. 
Pollution  factors  and  and  effectivity  indices  in  the  shrinking  mesh-cell 
Square  domain  shown  in  Fig.  5.2,  smooth  solution  u(X|,Z3)  =  4 sin  tX}  sinh A  = 
1(P,  p  =  2.  Note  that  the  pollution  factor  ocnverges  to  100%,  the  pollution  factor 
4^  diverges  to  infinity  and  the  effectivity  index  converges  to  zero  asymptotically; 
however  the  effectivity  index  converges  to  one  as  the  mesh-size  tends  to  zero. 


Pollution  factors  and  effectivity  indices  in  uff 


Exact  solution  u(xj,X3)  =  >4 sin xxj sinh xxj;  quadratic  elements 


Shrinking  mesh-cell  centered  at  (0.5,  0.5) 


Uniform  mesh  refinements  in  subdomain  only 


Mesh-size  in  the  rest  of  the  domain  =  0.125 


96.46 

99.76 

99.98 

365.79 

1463.63 

5851.72 

"o 


II 

■  in’".  Hi.* 

Effectivity  Index  1 

1.013 

0.998 

1.010 

1.000 

1.010 

0-999 

aiiB 


Figure  5<7.  Pollution-effect  for  meshes  refined  locnlly  in  nn  interior  subdomain.  The 
convex  domain  meshed  with  a  coarse  uniform  grid  with  mesh-size 

8 

The  mesh  has  been  divided  uniformly  fou'  times  in  an  interior  subdomain 


Qq  (the  mesh-size  in  the  subdomain  is  h  —  . 

\  128/ 
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Table  5.5a.  Pollution-effect  for  meshes  refined  locally  in  an  interior  subdomain.  Convex  domain  shown 
4 

in  Fig.  5.7,  a  =  p  =  2  (2a  >  p).  The  Vj  and  components  of  the  error  were  determined  from  the 
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5.5b.  Pollution-effect  for  meshes  refined  locally  in  an  interior  subdomain. 
Pollution  factors  and  and  effectivity  indices  in  the  shrinking  mesh-cell 
Ck>nvex  domain  shown  in  Fig.  5.7,  oi  =  ^,  p  =  2  (2a  >  p).  Note  that  the  pollution 
factor  converges  to  100%,  the  pollution  factor  diverges  to  infinity,  and  the 
effectivity  index  converges  to  zero.  However,  the  effectivity  index  converges  to 
one. 


Pollution  factors  and  effectivity  indices 


Shrinking  mesh-cell  centered  at  (0.5,  0.5) 

Exact  solution  u(r,0)  =  sin(y);  quadratic  elements 


Uniform  mesh-refinement  in  the  subdomain  only 


Mesh-size  in  the  rest  of  the  domain  h  ^  0.125 
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Dil'ccLionnI  Dorival.ivo  of  Iho  Rrror 


Study  of  the  Pollution  Errof^. 

Exact  SoluUon:  Un(r.e)  =  Sin(20/3) 

Quadratic  Elements;  Uniform  Refinements  in 
Mesh  Size  in  the  Subdomain  =  0.03125 
Mesh  Size  in  Rest  of  the  Domain  =  0.125 
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Figure  5.8.  Pollution-effect  for  meshes  refined  locally  in  an  interior  subdomain.  Quad¬ 
ratic  elements  (p  =  2),  a  =  Directional  derivatives  of  the  Vj  -  and 
Vj  -components  of  the  error  and  the  ZZ  error  estimator  along  the  s  axis  as 
shown  in  Fig.  3.1b.  (a)  Uniform  mesh-size  in  the  subdomain  h  «  (b) 

1  ^ '/}!/  I 

Uniform  mesh-size  in  the  subdomidn  h  - .  Note  that  — —  Z—L.\ 

Qy  128  ds  ^  ds\ 

and  asymptotically,  in  the  interior  of  the  subdomain. 
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Direclional  Derivative  of  the  Error 


5.2.S  Locally  refined  meshes  of  cubic  elements 


5.2.S.I.  Smooth  solution 

We  ocmsidered  the  same  boundary-value  problem  as  in  Section  5.2.1. 1  with 
the  smooth  soluti<m  tt(xi,X3)  s  Asinxxisinhrxj  and  employed  meshes  of  cubic 
elements  {p  —  2)  with  various  mesh-sises  in  the  subdomain.  In  Tables  5.6a,  5.6b  we 

r  the  p<dlution 
<c^.  We  note 

that,  as  the  mesh-size  io  the  subdomain  tends  to  zero,  the  pollution  factor 

converges  to  100%,  ^  diverges  to  infinity  and  the  effectivity  index  converges 
to  zero. 


give  the  values  of  |||V,  ||£^,  |||Va  |||^  and  j •  V, ,  the  values  o 
factors  and  and  the  values  of  the  effectivity  index  and 


5.2.S.2.  Singular  solution 

We  also  solved  the  mixed  boundary  value  problem  for  the  Li4>lacian  with 

boundary  conditions  consistent  with  the  exact  solution  u(r,6)  s  rtsin(7^)  and 

5 

employed  meshes  of  cubic  elements  (p  ~  3)  with  various  mesh-sizes  in  the  subdo¬ 
main  shown  in  Fig.  5.9.  In  Fig.  5.10a  (resp.  Fig.  5.10b)  we  show  the  gri^h  of  the 
directional  derivative  of  the  Vj  -  and  Vj  -components  of  the  error  in  the  subdomain, 

along  the  s-axis  shown  in  Fig.  3.1c,  where  the  mesh-size  in  the  subdomain  ^  ^  ~  ^ 

(resp.  h  ss  Once  again  we  note  that  as  the  mesh-size  in  the  subdomain  tends 

to  zero  the  ^-component  of  the  error  remains  constant  while  the  Vi  -component 
converges  to  zero.  Also  note  that  the  error  estimator  follows  the  -component 
of  the  error  and  that  the  -component  is  the  dominant  part  of  the  error  in  the 
interior  of  the  subdomain. 


5.3  Summary  of  the  results  for  the  pollution-effect  for 
locally  refined  meshes 


In  summary  we  observe  that  if  the  mesh  is  refined  locally  in  an  interior 
subdomain  only 


(i).  As  the  mesh  size  in  the  subdomain  tends  to  zero  the  Vj  -component  of  the  er¬ 
ror  converges  to  zero  in  the  interior  of  the  subdomain  while  the  -component 
remains  practically  constant. 


(ii).  The  value  of  relative  magnitude  of 


lliv;  iiu 


which  measures  the  poUution- 


effect  in  the  interior  of  the  subdomain  for  a  given  grid,  depends  on  the  mesh- 
size  in  the  subdomain  (compared  to  the  mesh-size  of  the  outside  mesh),  the 
smoothness  of  the  solution  and  the  degree  p  of  the  elements. 
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Tible  5.6b.  Pollution-effect  for  meshes  refined  locally  in  an  interior  subdomain. 
Pollution  factors  and  and  effectivity  indices  in  the  shrinking  mesh-cell 
Square  domain  shown  in  Fig.  5.2,  smooth  solution  u(zi,X3)  —  iisintrzjsinhs’Z],  A  — 
1(P,  p  s  3.  Note  that  the  pollution  factor  converges  to  100%,  the  pollution  fsctor 
diverges  to  infinity  and  the  effectivity  index  converges  to  zero  asymptotically; 
however  the  effectivity  index  k\  converges  to  one  as  the  mesh-size  tends  to  zero. 


Pollution  factors  and  effectivity  indices  in  uq 


Exact  solution  u(zi,Z])  =  Asinxzjsinhxz);  cubic  elements 


Shrinking  mesh-cell  u/q  centered  at  (0.5,  0.5) 


Uniform  mesh  refinements  in  subdomain  only 


Mesh-size  in  the  rest  of  the  domain  =  0.125 
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Figure  5.9.  Pollution-effect  Ibr  methei  refined  locally  in  an  interior  subdomain.  The  tri¬ 
angular  meshed  with  a  coarse  uniform  grid  with  mesh-aae  ^  —  g* 

The  mifh  has  been  divided  uniformly  four  times  in  an  interior  subdomain 
flo  (the  mesh-sise  in  the  subdomain  is  h  =  — ^ . 
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Directional  Derivative  of  the  Error 


l.OE-6 


-l.OE-6 


Studv  of  the  Pollution  Error^, 

Exact  Solution:  U|^r.6)  =  Sinf6e/5) 

Cubic  Elementa;  Uniform  Mesh  Refinements  in  Subdomain  only 
Mesh  Size  in  the  Subdomain  =  0.03125 
Mesh  Size  in  Reat  of  the  Domain  =  0.125 


_ ^ 


-3.0E-6  - 

; •»**•  Vi  Error  From  Overkill  (p  =  5) 

■  Vj  Error  From  Overkill  (p  =  5) 

Vi  Error  From  Smoothing  fp  =  6) 
Vt  Error  From  Smoothing  (p  -  6) 
ZZ  Error  Estimator 

-5.0E-6  4 'TT I  I’l'i’rri  i  in  i  ri'rri  f'ri'TrtTi”*  ih'T'H't  p-rt  ri  iTrn 
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Arc  Length  Coordinates  Along  S-Axis 


Figure  5.10.  Pollution-effect  for  grids  refined  locally  in  an  interior  subdomain.  Cubic 

elements  (p  =  3),  a  s  -.  Directional  derivatives  of  the  V^-  and  V^- 

o 

components  of  the  error  and  the  ZZ  error  estimator  along  the  a-axis  shown 
in  Fig.  3.1c.  (a)  Uniform  mesh-sise  in  the  subdomain  h  =  (b)  Uniform 

mesh-size  in  the  subdomain  h  =  Note  that  1^^  > 

12o  I  as  as 

dV 

(o’**  —  V«v)  •  a  »  asymptotically,  in  the  interior  of  the  subdomain. 

as 


DtrecliQnal  Derivative  of  the  Error 


Figure  5.10.  (continued) 
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6  Control  of  the  pollution-effect  in  a  subdomain 

Solutions  with  multiple  singularities  occur  cmnmonly  in  typical  2-D  geometrical 
idealizaticms  of  problems  in  solid-mechanics  and  heat  transfer.  Globally  adi4>tive 
methods  employ  nearly-equilibrated  grids  in  order  to  achieve  prescribed  uniform 
accuracies  over  the  entire  domain.  In  complex  engineering  problems  one  is  in¬ 
terested  to  control  the  error  and  to  compute  with  high-accuracy  only  in  certain 
subdomains  of  interest  (e.g.  discussion  about  global/local  analysis  in  the  introduc¬ 
tion).  Here  we  give  model  computations  which  demonstrate  that  it  is  possible  to 
guarantee  the  reliability  of  error  estimation  (and  to  control  the  error)  in  any  inte¬ 
rior  subdomain  Oq,  by  using  meshes  which  are  less  refined  outside  the  subdomain 
than  the  globally  adaptive  grids. 


6.1  Effect  of  the  refinement  outside  a  subdomain  on  the 
effectivity  index  in  the  interior  of  the  subdomain 

We  considered  the  mixed  boundary-value  problem  for  the  Laplacian  over  the 
Xr-shaped  domain  and  we  applied  boundary-conditions  consistent  with  the  exact 

solution  u(r,d)  :=  ri  sin(— d)  and  we  considered  a  subdomain  Qq  in  the  interior  of 

the  mesh,  far  from  the  singular  point.  We  refined  the  mesh  inside  the  subdomain 

^  III  III 

tiniformly,  to  a  desired  mesh-size  A,  and  computed  the  relative  error  Cn©  =  f||u  |j| 

using  a  pollution-free  finite  element  solution  over  the  subdomain.  (The  poUution- 
fne  solution  was  computed  by  solving  a  Neumann  problem  over  fig  using  data 
consistent  with  the  exact  solution.)  The  grid  outside  the  subdomain  fig  was  then 
generated  using  the  adaptive  refinement  algorithm  given  in  Section  4.5  with  a 
tolerance  S  =  where  A"  >  1  denotes  the  weight- factor.  In  the  following  we 

refer  to  meshes  of  this  type  as  weighted  equilibrated  meshes. 


6.1.1  Meshes  of  linear  elements 

In  Fig.  6.1  we  show  the  X-shaped  domain  meshed  by  a  grid  of  linear  elements 
(p  =  1).  The  grid  is  uniform  in  the  subdomain  Gg  :=  (0.5,0.75)*  with  mesh-size 

A  =  The  grid  outside  the  subdomain  was  generated  using  adaptive  refinem^t 

for  Coo  =  0.44  and  A  =  1.  In  Fig.  6.2  we  show  the  graph  of  the  effectivity  index, 
Kjt,  in  the  eight-element  mesh-cell  Wg  at  the  center  of  the  subdomiun  fig,  plotted 
against  the  weight-factor  A  for  different  mesh-sizes  in  the  subdomain.  Note  that 
the  curve  is  practically  invariant  with  respect  to  the  mesh-size  in  the  subdomain. 
Also  note  that  the  effectivity  index  Kk  remains  close  to  one  up  to  a  value  of  A  = 
2.5.  In  Fig.  6.3  we  show  the  weighted  equilibrated  mesh  of  linear  elements  where 
the  grid  outside  the  subdomain  is  generated  using  Coo  =  0-^  ^  ~  2.5. 


Figure  6.1.  Control  of  the  pollution-effect  in  a  subdomain  by  employing  weighted- 

equilibrated  meshes.  Linear  elements,  exact  solution  tt(r,  0)  =  sin(f );  the 

mesh-size  in  the  subdomain  fig  is  A  =  The  mesh  outside  fig  is  generated 

9a 

using  adaptive  refinement  for  ^o,  -  0.44  and  JiT  =  1.  The  effectivity  index 
in  the  eight-element  mesh-cell  at  the  center  of  the  subdomain  Qg  is 
0.92. 
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Effectivity  Index 


Study  of  the  Pollution  Error 

Singular  Solution:  \Xo{r,Q)  =  r‘^'Sin(©/3) 

Linear  Elements;  Adaptively  Refined  Mesh 
Location  of  Subdomain:  (0.5,0.75)x(0.5.0.75) 

A  Effect  of  Adaptive  Refinements  Outside  the  Subdomain 
3  on  the  Effecuvity  Index  in  the  Interior  of  the  Subdomain 


0.20  i 

q  »»♦■♦■<)  Mesh  Size  in  the  Subdomain  =  0.0625 

3  — — Mesh  Size  in  the  Subdomain  =  0.03125 

q  « >  Mesh  Size  in  the  Subdomain  =  0.015625 

0.00  1  I' I  I  I  I  I  I  I  I  I  I  I  I  I  I  T-i 'I'T  I  ri'i  I  r-i— 1—1  I  I  I  I  r  ri 

1.00  3.00  5.00  7.00 

Weight  Factor  (  K  ) 


Figure  6.2.  Control  of  the  pollution-dFect  in  a  subdomain  by  employing  weighted- 
equiKbrated  meshes.  Graph  of  the  elfectivity  index  in  the  eight-element 
niesh-cell  u/q  at  the  center  of  the  subdomain  G0  plotted  versus  the  weight 
factor  K  for  different  mesh-sizes  in  the  subdomain  of  the  shown  in 
Pig.  6.1a.  Note  that  the  curve  is  invariant  with  respect  to  the  mesh-sue  in 

(2q. 


134 


Figure  6.3.  Control  of  the  pollution-effect  in  a  subdomain  by  employing  weighted- 
equilibrated  meshes.  Linear  elements,  exact  solution  tt(r,^)  =  r^sin(|); 
the  mesh-size  in  the  subdomain  Qq  a  h  =  The  mesh  outside  the 
subdomain  is  generated  using  adaptive  refinement  for  Cn,  =  0*^  ^  — 

2.5. 
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In  Figs.  6.4&,  6.4b  and  6.4c  we  show  the  directional  derivatives  of  the  V|  -  and 
V2  'Components  of  the  error  along  the  s-axis  shown  in  Fig.  6.1a  for  K  =  1,  K  — 
2.5  and  if  =  5,  respectively.  Note  that  for  K  =  1,  the  directional  derivative  of 
the  V2  -component  of  the  error  is  close  to  zero  in  the  interior  of  the  subdomain. 
FHirther,  note  that  as  the  value  of  K  is  increased  i.e.  as  the  amount  of  refinement 
outside  the  subdomain  is  reduced,  the  directional  derivative  of  the  -component 
of  the  error  becomes  larger  in  the  intmor  of  the  subdomain  and  for  iif  =  5  the 
directional  derivative  of  the  -component  is  larger  than  the  directional  derivative 
of  the  Vi  -component  in  the  interior  of  the  subdomain.  Also  note  that  the  error 
estimator  practically  coincides  with  the  -component  of  the  error  in  the  interior 
of  the  subdomain. 

6.1.2  Meshes  of  quadratic  and  cubic  elements 

In  Fig.  6.5  we  show  the  Xrshaped  domain  meshed  by  a  grid  for  quadratic 
elements  (p  =  2).  The  grid  is  uniform  in  the  subdomain  fig  =  (0.5,0.75)^  with 

mesh-size  ^  ~  outside  the  subdomain  was  generated  using  adaptive 

refinement  for  ^  =  0.07  and  =  1.  In  Fig.  6.6  we  show  the  mesh  where  the  grid 
outside  the  subdomain  fig  was  generated  using  Cq,  =  0.01  and  K  =  1. 

In  Fig.  6.7  we  give  the  graph  of  the  effectivity  index,  ic^,  in  the  mesh-cell 
uf0  at  the  center  of  the  subdomain  fig  :=  (0.5,0.75)^  for  p  =  1,  2  and  3.  Note 
that  for  quadratic  elements  the  effectivity  index  is  close  to  one  up  to  a  value  of 
K  Z  and  for  cubic  elements  the  value  of  the  effectivity  index  is  close  to  one 
up  to  a  value  of  K  —  3.5.  In  Fig.  6.8  we  show  the  weighted-equilibrated  mesh 
of  quadratic  elements  which  is  generated  using  Cn,  =  0.07  and  K  =  Z  outside 
the  subdomain  and  in  Fig.  6.9  we  show  the  weighted-equilibrated  mesh  of  cubic 
elements  generated  using  adaptive  refinement  for  Coo  =  0.02  and  K  =  3.5  outside 
the  subdomain.  Hence  we  observe  that  as  the  polynomial  degree  p  of  the  elements 
is  increased,  the  pollution-effect  in  the  interior  of  the  subdomain  is  controlled  for 
a  bigger  range  of  values  of  the  weight-factor  K. 


6.2  Factors  aiOfecting  the  permbsible  range  of  the  weight- 
factor  K 

Here  we  examine  the  factors  which  affect  the  range  of  the  weight-factor  K  for 
which  the  pollution-effect  in  the  interior  of  the  subdomain  is  controlled.  We  will 
say  that  the  pollution-effect  in  an  interior  subdomain  is  controlled  if  >  0.85 

(where  u^g  denote  interior  mesh-cells  in  the  subdomain).  We  will  call  the  range 
of  K,  K  =  for  which  Kj,  >  0.85,  the  permissible  range  of  of  K.  For 

example,  in  the  numerical  examples  given  in  Section  6.1,  =  2.5,  3  and  3.5  for 

linear,  quadratic  and  cubic  elements,  respectively. 
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Directional  Derivative  of  the  Error 


0.006 


Location  of  the  Subdomain:  (6.5.0.75)x(0. 5.0.75) 
Weight  Factor  K  =  1;  Mesh  Size  in  the  Subdomain  = 


0.03125 


0.003 


-0.002 


-0.007 


■0.012  -TTT 

0.45 


•  Vi  Error  From  Overkill  (p  =  3) 
V*  Error  From  Overkill  (p  =  3) 
*■••■**22  Error  Estimator 


•  W  A • w 

Arc  Length  Coordinates  Along  S-Axis 


Figure  6.4.  Control  of  the  pollution-effect  in  a  subdomain  by  employing  weighted- 

equilibrated  meshes.  Linear  dements,  exact  solution  sin(|); 

the  mesh-ase  in  the  subdomain  Rq  is  A  =  Directional  derivatiyes  of 

the  Vi  -  and  -components  of  the  error  and  the  ZZ  estimator  along  the  a 

axis  for  adaptive  meshes  generated  using  (a)  =  1  (b)  K  —  2JS  {c)  K  — 

\8V  idF  I 

5.  Note  that  fox  K  =  I  and  K  =  2.5,  <.  but  for  JT  =  5, 

I  vs  I  vs  I 

^^1  in  the  interior  of  the  subdomain. 

as  \  as  \ 
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Directional  Derivative  of  the  Error 


0.008 


Study  of  the  Pollution  Errox^. 

Exact  Solution:  Uat(r.6)  =  Sin(6/3) 

Linear  Elements;  Adaptively  Refinea  Mesh 
Location  of  the  Subdomain:  (0.S.0.7S)x(0.5,0.75) 
Weight  Factor  K  =  2.5;  Mesh  Size  in  Uie  Subdomaii 


0.003 


-0.002  H 


-0.007 

H  \m 

Error  From  Overkill  (p  =  3) 
V,  Error  From  Overkill  (p  =  3) 
•••«•*«  22  Error  Estimator 
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Figure  6.4.  (continued) 


=  0.03125 
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Directional  Derivative  of  the  Error 


o.ooe 


Study  of  the  Pollution  Errox^, 

Exact  Solution:  u^r,6)  =  Sin(0/3) 

^  Linear  Elements;  Adaptively  Refined  Mesh 
Location  of  the  Subdomain:  (0.5.0. 75)x(0.5,0.75) 
Weight  Factor  K  s  5;  Mesh  Size  in  the  Subdomal 


0.003  i 


-0.002 


-0.007 


-0.012 


0.45  0.75  1.05 

Arc  Length  Coordinates  Along  S-Axis 
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Figure  6.4.  (continued) 


s  0.03125 
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Figure  6.5.  Control  of  the  pollution-effect  in  a  subdomain  by  employing  weighted- 
equilibrated  meshes.  Quadratic  elements,  exact  solution  =  ri  8m(|); 

the  mesh-size  in  the  subdomain  {lffi8h=  — .  The  mesh  outside  Oo  w  gen¬ 
erated  iiMtig  adaptive  refinement  for  Cn,  =  0*®^  K  —  1.  The  effectivity 
index  in  the  eight-element  mesh-cell  wj  at  the  center  of  the  subdomain 
CIq  is  0.97. 


Figure  6.6>  Control  of  the  pollution-effect  in  n  subdoniain  by  employing  weighted- 
equilibrated  meshes.  Cubic  dements,  exact  solution  tt(r,d)  =  sin(|);  the 

mesh-size  in  the  subdomain  (Ig  is  h  =  The  mesh  outside  Qg  is  generated 
using  adaptive  refinement  for  -  0.01  and  K  =  1.  The  effectivity  index 
in  the  dght-dement  mesh-cdl  at  the  center  of  the  subdomain  fig  is 
0.97. 


Effectivity  Index 


Fignxe  6.7.  Control  of  the  pollution-effect  in  a  lubdomain  by  employing  weighted- 
equilibrated  meihes.  Exact  lolution  ii(r,6)  =  ri  iin(|);  the  meih-nie  in 
the  lubdomain  Og  it  h  —  Compariion  of  the  effectivity  index  vi.  K 
curvet  for  p  =  1, 2, 3.  Note  that  at  the  polynomial  degree  p  of  the  elementi 
it  increated,  it  it  pottible  to  control  the  pdUntion-effect  in  the  tubdomain 
for  larger  valuei  of  the  weight  &ctor  K. 
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Figure  6.8.  Control  of  the  pollution-effect  in  a  subdomain  by  employing  weighted- 

equilibrated  meshes.  Quadratic  elements,  exact  solution  it(r,6)  =  ri  sin(|); 

the  mesh-size  in  the  subdomain  Qg  is  h  =  The  mesh  outside  Hg  is  gen- 

92 

erated  using  adaptive  refinement  for  Cn«  =  K  =  Z.  The  effectivity 

index  ic^  in  the  eight-element  mesh-cdll  at  the  center  of  the  subdomain 
^  Qg  is  0.93. 
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Figure  6.9.  Control  of  the  pollution-effect  in  a  subdomain  by  employing  weighted- 
equilibrated  meshes.  Cubic  dements,  exact  solution  u(r,ff)  =  sin(|);  the 

mesh-size  in  the  subdomain  (Ig  is  h  =  The  mesh  outside  Qg  is  generated 
using  adaptive  refinement  for  =  0.01  and  JC  =  3.5.  The  effectivity  index 
in  the  dght-dement  mesh-cdl  at  the  center  of  the  subdomain  fig  is 
0.93. 


144 


a.  Effect  of  the  distance  of  the  subdomain  from  the  singularity 

We  considered  the  Xrshaped  domain  with  a  uniform-mesh  subdomain  Hq 
(0.25,  0.5)’  closer  to  the  singular  point  than  the  subdomain  in  the  examples  given 
in  Section  6.1.  In  Fig.  6.10  we  show  the  mesh  for  linear  elements  (p  =  1)  for  which 

the  mesh-size  in  the  subdomain  ia  h  =  —  and  the  mesh  outside  the  subdomain 
was  generated  using  adaptive  refinement  for  Co,  =  0.25  and  /T  =  1.  In  Fig.  6.11 
we  show  the  graph  of  the  effectivity  index,  in  the  mesh-cell  Uq  at  the  center 
of  the  subdomain,  plotted  against  the  weight-factor  K.  Note  that  the  effectivity 
index  remains  close  to  one  up  to  a  value  of  /if  =  2  and  then  drops  to  small 
values  faster  when  compared  with  the  graph  of  Fig.  6.7.  In  Fig.  6.12  we  show  the 
weighted-equilibrated  mesh  of  linear  elements  where  the  grid  outside  the  subdomain 
is  generated  using  C  —  0.25  and  K  =  2.  Hence,  we  note  that  as  the  distance  of 
the  subdomain  Og  singular  point  increases,  the  permissible  range  of  the 

weight-factor  K  increases  in  size  and  the  pollution-effect  in  the  subdomain  can  be 
controlled  by  employing  coarser  meshes  outside  the  subdomain. 

b.  Effect  of  the  strength  of  the  singtUarity:  Value  of  the  exponent,  a 

We  considered  singularities  with  tvK>  different  exponents  ^  and  a  = 
for  quadratic  elements  and  plotted  the  graph  of  the  effectivity  index  for  tne 
mesh-cell  at  the  center  of  the  subdomain  (Ig  0.75)’  vs.  the  weight- 

factor  K.  In  Fig.  6.13  we  show  the  comparison  of  the  effectivity  index  vs  K  curves 

for  a  SB  -  and  ~  r.  From  the  graph,  we  note  that  =  3  for  a  =  while 
2  2 

K _ =  3.5  for  ^  ~  Hence  we  note  that  for  a  =s  -  the  permissible  range  of  the 

weight-factor  K  increases  in  size  and  the  values  of  the  effectivity  index  are  higher 
than  the  values  of  the  effectivity  index  for  or  s=  for  all  values  of  the  weight-factor 
K. 


c.  Effect  of  the  strength  of  the  singularity:  Value  of  the  stress  intensity  factor,  A 

We  considered  the  exact  solution  u(r,^)  =  Ar^  sin(^<>)  for  A  =  0.02,  A  =  0.1, 
A  =  0.5,  A  =  1,  A  =  10  and  A  =  50  and  linear  elements.  We  plotted  the  graph 
of  the  effectivity  index  for  the  mesh  cell  Uq  at  the  center  of  the  subdomain 
fig  :=  (0.5,  0.75)  X  (0,5,  0.75)  vs.  the  weight-factor  K.  In  Fig.  6.14  we  show  the 
comparison  of  the  graphs  for  A  :=  0.02,  0.1,  0.5,  1,  10  and  50.  Note  that  as 
the  value  of  stress-intensity  factor  A  is  increased,  the  value  of  the  effectivity  index 
drops  for  the  same  value  of  K.  Hence  we  observe  that  for  higher  values  of  the  stress 
intensity  factor  A,  the  permissible  range  of  the  weight-factor  K  decreases  and  more 
refinement  is  needed  outside  the  subdomain  in  order  to  control  the  pollution-effect 
in  the  interior  of  the  subdomain. 
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Figure  6.10.  Control  of  the  pollution-effect  in  a  subdomain  by  employing  weighted- 
equilibrated  meshes.  Effect  of  the  distance  of  the  subdomain  from  the 
singularity  on  the  amount  of  refinement  needed  outside  the  subdomain  in 
order  to  control  the  pollution-effect  in  Og*  Linear  elements,  exact  solution 
sin(f),  location  of  the  subdomain  Qg  :  (0.25,0.5)^,  mesh-sise  in 
the  subdomain  h=  The  mesh  outride  ftg  is  generated  using  adaptive 
refinement  for  (q,  =  0.25  and  =  1. 


Effectivity  Index 


Figure  6.11.  Control  of  the  pollution-effect  in  n  subdomain  by  employing  weighted- 
equilibrated  meshes.  Effect  of  the  distance  of  the  subdomain  (Iq  £rom  the 
singularity  on  the  amount  of  refinement  needed  outside  the  subdomain  in 
order  to  control  the  pollution-effect  in  fig.  Linear  elements,  exact  solution 
=  ri  sin(|),  location  of  the  subdomain  (Ig  :  (0.25,0.5)’,  mesh-sise 
in  the  subdomain  h  ~  The  mesh  outside  the  subdomain  is  generated 
using  adaptive  refinement  for  =  0.25  and  K  =  1.  Note  that  as  the 
distance  of  the  subdomain  from  the  nngnlar  point  increases,  it  is  posnUe 
to  control  the  pollution-effect  in  (2g  for  larger  values  of  the  weight-factor 
K. 
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Figtire  6.12.  Control  of  the  pollution-effect  in  a  subdomain  by  employing  weighted- 
equilibrated  meshes.  Effect  of  the  distance  of  the  subdomain  from  the 
singularity  on  the  amount  of  refinement  needed  outside  the  subdomain  in 
order  to  control  the  pollution-effect  in  the  interior  of  the  subdomain.  Linear 
elements,  exact  solution  v(r,d)  =  ri  8in(|),  location  of  the  subdomain  Qg  * 
(0.25,0.5)’,  mesh-size  in  the  subdomain  ^  The  mesh  outside  the 

subdomain  is  generated  using  adaptiye  refinement  for  Cn,  ==  0*25  and  K  — 

2. 
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Effectivity  Index 


Study  of  the  Pollution  Error 

Singular  Solution;  UEx(r.O)  =  r‘‘'“'Sin(0/3) 


Figure  6.13.  Control  of  the  pollution-effect  in  a  subdomain  by  employing  weighted- 
equilibrated  meshes.  Effect  of  the  exponent  a  on  the  amount  of  refinement 
needed  outside  outside  the  subdomain  in  order  to  control  the  pollution- 

effect  in  the  interior  of  the  subdomain.  Quadratic  elements,  exact  solution 

X  9  1 

«(r,5)  =  risin  mesh-size  in  the  subdomain  flo  :=  (0.5,0.75)^,  h  =  — . 

Note  that  as  the  value  of  the  exponent  increases  (i.e.  the  solution  becomes 

smoother),  it  is  possible  to  control  the  pollution-effect  in  Qq  for  larger  values 

of  the  wdght-factor  K. 


149 


Effectivity  Index 


Figure  6.14.  Control  of  the  pollution-effect  in  a  subdomain  by  employing  weighted- 

equilibrated  meshes.  Effect  of  the  stress-intensity  factor  A  on  the  amount  of 

refinement  needed  outside  the  subdommn  in  order  to  control  the  pollution- 

effect  in  the  interior  of  the  subdomain.  Linear  elements,  exact  solution 
id 

«(r,d)  =  rJ  sin  location  of  the  subdomain  SIq  :=  (0.5,0.75)’,  mesh-size 

O  ^ 

in  the  subdomain  h  =  — .  Note  that  as  A  is  decreases,  it  is  possible  to 
control  the  pollution-effect  in  Qg  for  larger  values  of  the  weight-factor  K. 
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Effectivity  Index 


Figure  6.15.  Control  of  the  pollution-effect  in  a  subdomain  by  employing  weighted- 
equilibrated  meshes.  Smooth  solution  ti(xj,Z2)  =  sin  ira;^  sinh  tZ]!  mut 
square  Hnnutin  shown  in  Fig.  5.2.  Comparison  of  the  effectivity  index  vs.  K 
curves  for  p  =  1,  2,  3.  Note  that  as  p  increases,  it  is  possible  to  control  the 
pollution-effect  in  (Ig  piacticatty  for  all  values  of  the  weight-factor  K. 
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d.  Effect  of  the  smoothness  of  the  solution 

We  also  considered  the  smooth  solution  u(xi,Z3)  s=  sin irZ}  sinh irx,  in  the 
square-domain  =  (0,1)^  and  plotted  the  graph  of  the  effectivity  index  ic^, 

for  the  mesh  cdl  at  the  center  of  the  subdomain  fig  (0.375,0.625)’  as  shown 
in  Fig.  5.2,  vs.  the  weight  factor  K.  In  Fig.  6.15  we  show  the  graphs  for  p  =  1,  2 
and  3.  Note  that  the  values  of  the  effectivity  index  are  greater  than  0.8  for  the  all 
values  of  K. 

6.3  Summary  of  the  results  for  the  control  of  the  pollution- 
effect  in  a  subdomain 

In  summary,  we  observe  that 

(i) .  It  is  possible  to  control  the  pollution-effect  in  an  interior  subdomain  by  using 

weighted-equilibrated  grids  which  are  less  refined  outside  the  subdomain  than 
globally  nearly  equilibrated  grids. 

(ii) .  The  amount  of  refinement  outside  the  subdomain  needed  to  control  the 

pollution-effect  in  the  subdomain  depends  on  the  strength  of  the  singulari¬ 
ties,  the  distance  of  the  subdomain  fig  from  the  singular  points,  the  mesh-size 
inside  the  subdomain  and  the  polynomial  degree  p  of  the  elements. 
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7  Summary 

In  thin  work  we  presented  a  numerical  study  of  the  pollution*error  in  the  k- 
version  of  the  finite^ement  method.  We  studied  the  effect  of  the  pollution  error 
on  the  local  quality  of  error  estimation  for  various  classes  of  meshes  and  solutions 
used  in  practical  engineering  computations.  In  particular,  we  considered  harmonic 
solutions  with  algebraic  singularities  of  the  type  and  studied  the  pcdlution-effect 
for  uniform  meshes,  globally  adaptive  (nearly  equilibrated)  meshes,  locally  refined 
meshes  and  weighted-equilibrated  meshes. 

The  major  conclu8i<m8  of  this  study  are: 

1.  When  the  mesh  is  globally  adaptive  (nearly  equilibrated  in  the  energy-norm) 
the  pollution-effect  is  negliffble  and,  for  all  practical  purposes,  we  can  guar¬ 
antee  the  reliability  of  error  estimation  for  all  interior  mesh-cells. 

2.  If  the  mesh  is  not  nearly  equilibrated  in  the  energy-norm,  the  pollution-effect 
may  be  significant. 

3.  For  uniform  meshes,  the  poUution-^ect  is  significant  asymptotically  when 
2a  <  p,  the  polynomial  degree  of  the  elements.  However  when  2a  >  p,  the 
pollution  is  negligible. 

4.  The  accuracy  of  the  finite  element  solution  in  a  subdomain  depends  upon  the 
relationship  between  the  mesh  inside  the  subdomain  and  the  mesh  outside 
the  subdomain  and  the  smoothness  of  the  exact  solution. 

5.  It  is  possible  to  control  the  pollution  effect  in  a  subdomain  by  employing 
weighted-equilibrated  grids  which  are  less  refined  outside  the  subdomain  than 
the  globally  nearly  equilibrated  grids. 

6.  The  reliability  of  local  error  estimation  in  an  interior  sub^^omain  can  be 
guaranteed  only  when  the  pollution-effect  in  the  subdomain  is  negligible.  If 
there  is  strong  pollution,  severe  underestimation  may  occur  and  in  general 
nothing  can  be  said  about  the  reliability  of  the  local  results  of  a-postenori 
error  estimation. 
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foose  in  foe  mathematical  community. 

To  provide  a  limited  consulting  service  in  all  areas  of  numerical  mathematics  to  foe 
University  as  a  whole,  and  also  to  government  agencies  and  industries  in  the  State  of 
Mai^and  and  foe  Washington  Metropolitan  area. 

To  wifo  foe  education  of  numerical  anatysts,  eq)cdally  at  foe  postdoctoral  level, 
in  conjunction  wifo  foe  Interdisc^linary  ^>plicd  Mathematics  Program  and  foe 
programs  of  foe  Mathematics  and  C^puter  S^noe  Departments.  This  includes  active 
collaboration  wifo  government  agencies  such  as  the  National  Institute  of  Standards  and 
Tedmology. 

To  be  an  international  center  of  study  and  researdi  for  foreign  students  in  numerical 
mathematics  who  are  supported  by  foreign  governments  or  ezdiaoge  agencies 
(Fulbright,  etc). 

Further  information  nuy  be  obtained  from  Professor  L  BabuSla,  Qiairman,  Laboratory  for 
Numerical  Analysis,  Institute  for  Physical  Sdenee  and  Tedmology,  University  of  Maryland,  CoUege 
Park,  Maryland  20742-2431. 


